{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Logistic Regression"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Why"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-02-25T01:47:41.319391Z",
     "start_time": "2019-02-25T01:47:40.945260Z"
    }
   },
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "import warnings\n",
    "\n",
    "warnings.simplefilter(\"ignore\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "描述一个事件的可能性，除了用概率以外，还可以用几率，胜率，赔率等等。"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-02-25T01:47:41.327409Z",
     "start_time": "2019-02-25T01:47:41.322549Z"
    }
   },
   "outputs": [],
   "source": [
    "p = np.linspace(0,1,100)\n",
    "odds = p/(1-p)\n",
    "logodds = np.log(odds)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-02-25T01:47:41.576818Z",
     "start_time": "2019-02-25T01:47:41.330941Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEKCAYAAAAIO8L1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAH7FJREFUeJzt3XuUXGWZ7/HvU1V979xIQgMJkkQCEiKXpAkwAiagDAeVgMfjXQHRuEYEjnhccmTNkjmzcMTxIM7oGcwIDKBjRI4Cig4ymJYThQjhJkkTCYGQQCCdkE66093VXVXP+WPv7lQ6laS7k9q7Lr/PWrVqX95d+3mSsB/2++6LuTsiIiLDJeIOQERESpMKhIiIFKQCISIiBalAiIhIQSoQIiJSkAqEiIgUpAIhIiIFqUCIiEhBKhAiIlJQKu4ADsaUKVN8xowZY9p2165dNDU1HdqAykA15l2NOUN15l2NOcPo8161atVWd596oHZlXSBmzJjBk08+OaZt29raWLhw4aENqAxUY97VmDNUZ97VmDOMPm8z2zCSdupiEhGRgopWIMzsdjPbYmbP5y07zMweNrMXw+9J4XIzs38ys3Vm9pyZzStWXCIiMjLFPIP4N+CCYcuuAx5x99nAI+E8wH8BZoefJcC/FDEuEREZgaKNQbj7o2Y2Y9jixcDCcPpOoA34arj8Lg+ePf64mU00syPdffNo9zswMMCmTZvo6+vbb7sJEybQ3t4+2p8vKfX19UyfPp2ampq4QxGRChT1IHVL3kH/DaAlnJ4GbMxrtylctleBMLMlBGcZtLS00NbWtsf65uZmWlpamDZtGma2z0Cy2SzJZHKMacTP3dmxYwfPPvss3d3dI96uu7t7rz+zSleNOUN15l2NOUPx8o7tKiZ3dzMb9duK3H0psBSgtbXVh4/ct7e3M3369P0WB4Curi7GjRs32t2XlHHjxtHd3U1ra+uIt6nGqzyqMWeozryrMWcoXt5RX8X0ppkdCRB+bwmXvwYcndduerhsTA5UHCpFteQpIvGIukA8AFwaTl8K3J+3/NPh1UxnADvGMv4gIlLpevozfPuhtTyzsbPo+yrmZa4/AR4DjjezTWZ2BfBN4L1m9iLwnnAe4NfAemAd8K/AF4oVV6m44YYb+Pa3v73X8ldeeYW5c+fGEJGIlIOdvRm+t3wda17fWfR9FfMqpo/tY9V5Bdo6cGWxYhERqRTpTBaAulTxO4B0J3WR3HzzzcydO5e5c+dyyy23AHDjjTdy3HHHcdZZZ7F27dqhtqtWreLkk0/m5JNP5vvf//7Q8tWrV7NgwQJOOeUUTjrpJF588cXI8xCR0pLO5ACoqyn+4busn8V0IH/3y9X7PA0b62Wuc44az9c/cOJ+26xatYo77riDlStX4u6cfvrpnH322SxbtoxnnnmGTCbDvHnzmD9/PgCXX3453/ve9zjnnHP4yle+MvQ7t956K9dccw2f+MQn6O/vJ5vNjjpeEaks6YGwQKSKf5m+ziCKYMWKFVxyySU0NTXR3NzMBz/4QR588EEuueQSGhsbGT9+PBdddBEAnZ2ddHZ2cs455wDwqU99auh3zjzzTL7xjW9w0003sWHDBhoaGmLJR0RKR5RdTBV9BrG//9Mvh/sgPv7xj3P66afz4IMPcuGFF/KDH/yAc889N+6wRCRGQ11MGoMoT2effTb33XcfPT097Nq1i1/84he8733v47777qO3t5euri5++ctfAjBx4kQmTpzIihUrAPjxj3889Dvr169n1qxZXH311SxevJjnnnsulnxEpHQMnUHUFL+LqaLPIOIyb948LrvsMhYsWADAZz/7WebPn89HPvIRTj75ZA4//HBOO+20ofZ33HEHn/nMZzAzzj///KHl99xzD3fffTc1NTUcccQRfO1rX4s8FxEpLf0RnkGoQBTJtddey7XXXrvHsuuvv57rr79+r7bz58/n2WefHZr/1re+BcB1113Hddddt1d7Eale6mISEZGChq5iiqCLSQVCRKSM6Ea5gxTcmF35qiVPEdlNXUwHob6+nm3btlX8wdPd2bZtG/X19XGHIiIR2l0gdBXTqE2fPp1NmzbR0dGx33Z9fX1lf3AdfKOciFSP9EAWM6hJFv9x/xVXIGpqapg5c+YB27W1tXHqqadGEJGIyKGTzuSoSyUieR9MxXUxiYhUsqBARPO6ZBUIEZEyks5kIxmgBhUIEZGykh7IRfKob1CBEBEpK+piEhGRgtTFJCIiBQ1exRQFFQgRkTKSHshRqwIhIiLDBV1MGoMQEZFh1MUkIiIFpTO5SB71DSoQIiJlJT2gq5hERKQAdTGJiEhBulFOREQKSmeyetSGiIjsKZtzBrKuLiYREdlTf4RvkwMVCBGRspHOZIFo3kcNKhAiImVj6H3UlTwGYWZfMrPVZva8mf3EzOrNbKaZrTSzdWb2UzOrjSM2EZFSlR6o8C4mM5sGXA20uvtcIAl8FLgJ+I67HwtsB66IOjYRkVJWLV1MKaDBzFJAI7AZOBe4N1x/J3BxTLGJiJSkoS6mSi0Q7v4a8G3gVYLCsANYBXS6eyZstgmYFnVsIiKlbPcYRDRdTKlI9pLHzCYBi4GZQCfwM+CCUWy/BFgC0NLSQltb25ji6O7uHvO25awa867GnKE68670nNu3BV1M7c8/h7++u0gUK+/ICwTwHuBld+8AMLOfA+8CJppZKjyLmA68Vmhjd18KLAVobW31hQsXjimItrY2xrptOavGvKsxZ6jOvCs9Z1+7BZ54gtNb53Hq2yYNLS9W3nGMQbwKnGFmjWZmwHnAGmA58KGwzaXA/THEJiJSsir+KiZ3X0kwGP0U8OcwhqXAV4FrzWwdMBm4LerYRERK2dBVTBHdBxFHFxPu/nXg68MWrwcWxBCOiEhZqPirmEREZGzSehaTiIgUkh6ItotJBUJEpEyoi0lERAoaLBC1SRUIERHJk85kqUslCO4QKD4VCBGRMpEeyFEbUfcSqECIiJSNdCYX2RVMoAIhIlI2BruYoqICISJSJtKZXGSXuIIKhIhI2UgPqItJREQKUBeTiIgUFAxSq0CIiMgwwRiEuphERGSY9IC6mEREpIB+dTGJiEghulFOREQKSmeyug9CRET2FtwHoQIhIiLDqItJRET2kss5/VmdQYiIyDD92fBtchqDEBGRfLtfN6ouJhERyZPOZIHo3kcNKhAiImUhPTB4BqECISIieYa6mPQsJhERyacuJhERKWj3ILUKhIiI5Nk9BqEuJhERyTPUxaT7IEREJJ+6mEREpKCqKRBmNtHM7jWzF8ys3czONLPDzOxhM3sx/J4UR2wiIqUoPTB4FVPlj0F8F/gPd38HcDLQDlwHPOLus4FHwnkREaFKziDMbAJwDnAbgLv3u3snsBi4M2x2J3Bx1LGJiJSqankW00ygA7jDzJ42sx+aWRPQ4u6bwzZvAC0xxCYiUpLiuIrJ3D2ynQGYWSvwOPAud19pZt8FdgJXufvEvHbb3X2vcQgzWwIsAWhpaZm/bNmyMcXR3d1Nc3PzmLYtZ9WYdzXmDNWZdyXn/IsX+7n/pQFu/+tGEmZ7rBtt3osWLVrl7q0HbOjukX6AI4BX8ubPBh4E1gJHhsuOBNYe6Lfmz5/vY7V8+fIxb1vOqjHvaszZvTrzruSc/+HX7T77a78uuG60eQNP+giO15F3Mbn7G8BGMzs+XHQesAZ4ALg0XHYpcH/UsYmIlKp0JhvpADVAKtK97XYV8GMzqwXWA5cTjIfcY2ZXABuAD8cUm4hIyUlncpGOP0BMBcLdnwEK9X+dF3UsIiLlID2Qi/QKJtCd1CIiZSGOLiYVCBGRMpDO5KhVgRARkeGCMQh1MYmIyDDpAXUxiYhIAelMTgVCRET2FhQIdTGJiMgw6Uw28vsgRrQ3M3u7mdWF0wvN7Gozm3ig7URE5NAI7oMowQIB/F8ga2bHAkuBo4F/L1pUIiKyh1LuYsq5ewa4BPhnd/8KwQP1REQkAv0lfKPcgJl9jOAher8Kl9UUJyQRERkujmcxjXRvlwNnAje6+8tmNhO4u3hhiYjIIHePpYtpRA/rc/c1wNV58y8DNxUrKBER2a0/G/37qOEABcLM/gzs85Vz7n7SIY9IRET2sPt91CVUIID3h99Xht+D3UqfZD+FQ0REDp30QFggIn4W034LhLtvADCz97r7qXmrvmpmTwHXFTM4EREJbpIDqEuW5iC1mdm78mb+ahTbiojIQRjqYirRN8pdAdxuZhMAA7YDnylaVCIiMmSoi6nExiAAcPdVwMlhgcDddxQ1KhERGbKzbwCAprpo3xJ9oKuYrt3HcgDc/eYixCQiInm2dqcBOHxcfaT7PVA5Ghd+Hw+cBjwQzn8A+FOxghIRkd06uoICMXVcXaT7PdBVTH8HYGaPAvPcvSucvwF4sOjRiYgIHV1pUgljYkO0Tzga6YhHC9CfN98fLhMRkSLr6EozpbmORMIi3e9IRzzuAv5kZr8I5y8B7ixOSCIikq+jOx159xKM/CqmG83sN8D7CO6gvszdny5qZCIiAgRnEC3jox2ghpG/Ue5qgjOGFFAL3GlmVxUzMBERCQRdTLWR73ekXUyfBc5w910AZnYT8Bjwz8UKTEREIJdztu3qj6WLacSP2gCyefPZcJmIiBTR9p5+sjlnanOJjkEAdwAr8wapLwZuK05IIiIyqKN78B6I6McgRjpIfbOZtQFnhYsu1yC1iEjxxXWTHIz8DAJ3fwp4qoixiIjIMHEWCD2yW0SkhFVlgTCzpJk9bWa/CudnmtlKM1tnZj81s+iv6RIRKTEdXWkaapI01Ub7NjmI9wziGqA9b/4m4DvufizB+yauiCUqEZESMngX9eBTtKMUS4Ews+kEd2X/MJw34Fzg3rDJnQRXSomIVLWtMT1mA8DcPfqdmt0L/APB48T/B3AZ8Hh49oCZHQ38xt3nFth2CbAEoKWlZf6yZcvGFEN3dzfNzc1j2racVWPe1ZgzVGfelZjz9St6OKIpwVWn7vsy19HmvWjRolXu3nqgdtG+nggws/cDW9x9lZktHO327r4UWArQ2trqCxeO+icAaGtrY6zblrNqzLsac4bqzLsSc+559LecMPNIFi585z7bFCvvyAsE8C7gIjO7EKgHxgPfBSaaWcrdM8B04LUYYhMRKRn9mRzbewaY2hz9TXIQwxiEu/9Pd5/u7jOAjwK/c/dPAMuBD4XNLgXujzo2EZFSsm1XfJe4QmndB/FV4FozWwdMRo/yEJEqF+c9EBBPF9MQd28D2sLp9cCCOOMRESklcReIUjqDEBGRPFu7VSBERKSAwTOIOF4WBCoQIiIlq6MrzYSGGupS0T9mA1QgRERKVkeMd1GDCoSISMmK613Ug1QgRERKVEdXOpY3yQ1SgRARKVEdXelY3kU9SAVCRKQE9fRn2NWf1RiEiIjsaWtXPxDfPRCgAiEiUpI6uvsAFQgRERlm0/ZeAFrGq0CIiEie9s1d1CSNWVPiewGSCoSISAlq37yT2YePozYV32FaBUJEpASt2byTE44cH2sMKhAiIiVmS1cfHV1p5hylAiEiInnaN3cBMEdnECIikm/N6zsBFQgRERmmffNOpk1sYEJjTaxxqECIiJSYUhigBhUIEZGS0tufZX1Hd+wD1KACISJSUta+2UXO4x9/ABUIEZGSMjhAfaLOIEREJF/75p2Mq0sxfVJD3KGoQIiIlJLBAWozizsUFQgRkVKRyzntm3eWxAA1qECIiJSMDW/10NOfLYkBalCBEBEpGYMD1KVwDwSoQIiIlIwV6zpoqk1y3BHxvQMinwqEiEgJyOWch9dsYeE7DqculYw7HEAFQkSkJDy9sZOt3WnOn9MSdyhDVCBERErAb9e8QSphLDz+8LhDGRJ5gTCzo81suZmtMbPVZnZNuPwwM3vYzF4MvydFHZuISFweXv0mZ759MhMa4n2Ca744ziAywJfdfQ5wBnClmc0BrgMecffZwCPhvIhIxVu3pZv1W3eVVPcSxFAg3H2zuz8VTncB7cA0YDFwZ9jsTuDiqGMTEYnDb9e8AcB7qr1A5DOzGcCpwEqgxd03h6veAErrT0pEpEh+u/pNTpo+gSMnxP/8pXzm7vHs2KwZ+D1wo7v/3Mw63X1i3vrt7r7XOISZLQGWALS0tMxftmzZmPbf3d1Nc3NpXGscpWrMuxpzhurMuxxz7uzL8d/bevng7BouenvtmH5jtHkvWrRolbu3HrChu0f+AWqAh4Br85atBY4Mp48E1h7od+bPn+9jtXz58jFvW86qMe9qzNm9OvMux5zveuwVP+arv/IXNu8c82+MNm/gSR/BsTqOq5gMuA1od/eb81Y9AFwaTl8K3B91bCIiUcrlnLv++ArvOGIcx7WU3plPHGMQ7wI+BZxrZs+EnwuBbwLvNbMXgfeE8yIiFet3L2zhxS3dfP7ds0ri8d7DpaLeobuvAPb1J3FelLGIiMTp1t+/xLSJDbz/pKPiDqUg3UktIhKDJ155iyc3bOdzZ8+kJlmah+LSjEpEpMLd2vYSkxpr+PBpR8cdyj6pQIiIRGztG1088sIWLvurmTTWRt7TP2IqECIiEbv54bU01CT59JnHxB3KfqlAiIhE6Dd/3sxDq9/ki+cey6Smsd0YFxUVCBGRiGzf1c/f3r+audPGs+ScWXGHc0Cl2/klIlJh/v5Xa+js6eeuzywo2SuX8pV+hCIiFWD5C1v4+dOv8YWFb2fOUePjDmdEVCBERIrsla27+PLPnuW4lmauPPfYuMMZMRUIEZEiemtXP5fd8SfcnaWfaqUulYw7pBHTGISISJH0DWT53F1P8vqOPn7yudOZMaUp7pBGRWcQIiJFkM5kuWbZ0zz16nZu+cgpzD/msLhDGjWdQYiIHGI7egf4/N1P8vj6t7jhA3O48J1Hxh3SmKhAiIgcQm/s6OOyO/7ESx3d3PKRU7j41GlxhzRmKhAiIofIYy9t40s/fYbudIY7LlvAWbOnxB3SQVGBEBE5SP2ZHN/5z79w6+9fYubkJm6/7LSyuddhf1QgREQOwjMbO/nb+57nz6/t4GMLjuZv3z+npJ/QOhqVkYWISMS2dPXxj/+xlp+t2sTUcXXc+sl5XDC3PAej90UFQkRkFLZ1p7ltxcvc9dgG0pksn3/3LK46dzbNdZV3OK28jEREimDjWz3c9dgr/OjxV+nLZLlw7pF8+fzjmDW1Oe7QikYFQkRkH7I559G/dPCjxzfwu7VbMGDxKdO4ctHbOfbwcXGHV3QqECIiedyd5zbt4P5nXueXz71OR1eaKc11fHHRsXxswds4amJD3CFGRgVCRKpefybHn15+i4fXvMF/tm/htc5eapMJFr1jKhefMo3zTmihNlV9TyZSgRCRqpPNOWvf6OKPL21lxbqtrFz/Fr0DWeprEpx17FSuOW82f33iEUxorIk71FipQIhIxevs6ee5TTt4ZmMnqzZs56kN2+lKZwCYNbWJ/9Y6nbNnT+WsY6fQUFs+j+MuNhUIEakYOXc2vtXD2je6aN+8kzXhZ8O2nqE2sw9v5gOnHMWCGYdx2szDmFZFYwqjpQIhImWnq2+ADdt62LCth/Ud3azfuouXOrpZu7mH9EPLh9rNmNzIiUeN58OtR3Pq0ROZO30C4+uru9toNFQgRKSkuDs7ezNs3tnL6529vNbZx+udvWx8q4eN23vZ9FYP23b177HNURPqmTW1mXOmpzh3/gkc19LM8UeMr8ib16KkPz0RicRANsf2Xf1s7e5na3eard1pOrqCz5auNG/u7GNLV5rNO3rpG8jtsW1N0jhqYgNHT2rkvXNaOGZyE8dMbuSYyY3MnNI09OyjtrY2Fi54WxzpVSQVCBEZlWzO6eobYGdvhh29A+zoHaCztz/47hmgs6ef7T3B9Paeft7aFXx29A4U/L36mgQt4+tpGVfPnKPGc947DueICfUcMaGeoyY2MG1iA1Oa60gmLOJMRQVCpAr0Z3L09mfpGcjQ05+lJ51lV3+GXekM3elg2eD04HdX3+7vrr6B8DtYtj/1NQkmNdYyoaGGw5pqOfGo8UxqrGVycy2Tm+uY3FTL1HF1TGmuY0pzLc11Kcx08C9FJVUgzOwC4LtAEvihu38z5pBEDrlczunP5khncvRncqQz2fA7/Axkh6b7wum+gexe030DwXRvON8bLuvpD+a3d/WQbXuI3v4smZyPOL6GmiTj6lM016dorksxrj7FlOYmxtXXML6+hvENKcbV1zChIfiMr08xsbGWiY3BfH2NLhOtFCVTIMwsCXwfeC+wCXjCzB5w9zXxRialyN3JedCvnc05mZwH39kcmZyTyTqZXDD96s4sz23qZCAbrM/mnIGcM5DJkcnlguWD31lnIJtjIPydgUwuaJvNBdPZHP3h7wTtgoP9QDY42A9+92ed/kyW/sH5wU+4zcGqr0lQX5OkIfzU1SRprA2mJzXW0FCbYse2NLPeNp2G2iRNtUkaalPhd5LG2hRNdUmaBr/rUsGnNqWuHBlSMgUCWACsc/f1AGa2DFgMlF2BcHfcg2uyc+H37vlg2eABLpvzoenB9YNts7ndbbPu5HJ7/sZe24bts0PTu39jcLvnX8+w/elN5HK72w1957XN7rF8WNv85bkcWQ+mB9dnc3tun8nlfXtwEB5cn78uGx7Qh2+TDQ/Q2bzfz+SCP6cR++MfDurvtCZp1CYTpJKJoemaVIJUwqhNJalNGjXJBLWpBI21KWpTwXTtYPtUgtpkktpUgrpwXd3QZ/fyuppgvi4VFIDB9fU1CeqG5hMj6pJpa2tj4cITDypvqW6lVCCmARvz5jcBpxdjR/c8sZFb/l8P9avahg7G+d97HrDDAz67D8L7O/g7jO7AFYfnnh3TZgmDZMJIWPBJJnZ/gnlIJRIkEpA0I5EwkmG7VHL39GD7hlQybAOpZGLPtuG2qWRwEB7cLpW0YD5cl0wE8/ntapJGMpEIlxtr29dwyknv3KtNTTJBKhEcwAfX1aYSQ783eHBPJkx95FKVSqlAjIiZLQGWALS0tNDW1jbq39j0ZoaWhhw1yT4Gz6bNIEEwk7Bgfmia8GMMtU+YYVi4XbDOCBomyG+3e12wDXtuM2z7/PZmNrQuYbtjGf6bQzEOa7Pn8qBdb08PzU2Ne7XL/52E7f0ZjGf0PPwUkQP7GTc9vqmP5JvtODAQfvqKG1FJ6O7uHtN/H+WsGnOG4uVdSgXiNeDovPnp4bI9uPtSYClAa2urL1y4cNQ7WgjMa2tjLNuWu7YqzLsac4bqzLsac4bi5V1Kz699AphtZjPNrBb4KPBAzDGJiFStkjmDcPeMmX0ReIjgMtfb3X11zGGJiFStkikQAO7+a+DXccchIiKl1cUkIiIlRAVCREQKUoEQEZGCVCBERKQgFQgRESnIvOSfC7FvZtYBbBjj5lOArYcwnHJRjXlXY85QnXlXY84w+ryPcfepB2pU1gXiYJjZk+7eGnccUavGvKsxZ6jOvKsxZyhe3upiEhGRglQgRESkoGouEEvjDiAm1Zh3NeYM1Zl3NeYMRcq7ascgRERk/6r5DEJERPaj4guEmV1gZmvNbJ2ZXVdgfZ2Z/TRcv9LMZkQf5aE1gpyvNbM1ZvacmT1iZsfEEeehdqC889r9VzNzMyv7q11GkrOZfTj8+15tZv8edYzFMIJ/428zs+Vm9nT47/zCOOI8lMzsdjPbYmbP72O9mdk/hX8mz5nZvIPeafD+5Mr8EDw2/CVgFlALPAvMGdbmC8Ct4fRHgZ/GHXcEOS8CGsPpvyn3nEead9huHPAo8DjQGnfcEfxdzwaeBiaF84fHHXdEeS8F/iacngO8EnfchyDvc4B5wPP7WH8h8BuCF0CeAaw82H1W+hnEAmCdu693935gGbB4WJvFwJ3h9L3AeVbeLyA+YM7uvtzde8LZxwne3lfuRvJ3DfD3wE1UxltHR5Lz54Dvu/t2AHffEnGMxTCSvB0YH05PAF6PML6icPdHgbf202QxcJcHHgcmmtmRB7PPSi8Q04CNefObwmUF27h7BtgBTI4kuuIYSc75riD4v45yd8C8w1Puo939wSgDK6KR/F0fBxxnZn8ws8fN7ILIoiuekeR9A/BJM9tE8I6Zq6IJLVaj/W//gErqhUESLTP7JNAKvDvuWIrNzBLAzcBlMYcStRRBN9NCgjPFR83sne7eGWtUxfcx4N/c/X+b2ZnA3WY2191zcQdWTir9DOI14Oi8+enhsoJtzCxFcDq6LZLoimMkOWNm7wGuBy5y93REsRXTgfIeB8wF2szsFYI+2gfKfKB6JH/Xm4AH3H3A3V8G/kJQMMrZSPK+ArgHwN0fA+oJnldUyUb03/5oVHqBeAKYbWYzzayWYBD6gWFtHgAuDac/BPzOwxGfMnXAnM3sVOAHBMWhEvqk4QB5u/sOd5/i7jPcfQbB2MtF7v5kPOEeEiP5930fwdkDZjaFoMtpfZRBFsFI8n4VOA/AzE4gKBAdkUYZvQeAT4dXM50B7HD3zQfzgxXdxeTuGTP7IvAQwZUPt7v7ajP7X8CT7v4AcBvB6ec6ggGgj8YX8cEbYc7/CDQDPwvH419194tiC/oQGGHeFWWEOT8EnG9ma4As8BV3L+cz5JHm/WXgX83sSwQD1peV+f/4YWY/ISj2U8Kxla8DNQDufivBWMuFwDqgB7j8oPdZ5n9mIiJSJJXexSQiImOkAiEiIgWpQIiISEEqECIiUpAKhIiIFKQCISIiBalAiIhIQSoQIoeQmc0wsxfM7Mdm1m5m95pZY9xxiYyFCoTIoXc88H/c/QRgJ8E7R0TKjgqEyKG30d3/EE7/CDgrzmBExkoFQuTQG/78Gj3PRsqSCoTIofe28B0EAB8HVsQZjMhYqUCIHHprgSvNrB2YBPxLzPGIjElFP+5bJCYZd/9k3EGIHCydQYiISEF6H4SIiBSkMwgRESlIBUJERApSgRARkYJUIEREpCAVCBERKUgFQkRECvr/8hB0X7iuf4sAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(p, odds, label=\"odds\")\n",
    "plt.xlabel(\"p\")\n",
    "plt.ylabel(\"odds\")\n",
    "plt.legend()\n",
    "plt.grid()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-02-25T01:47:42.218499Z",
     "start_time": "2019-02-25T01:47:41.579902Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEWCAYAAABxMXBSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xl8VeW59//PRUjCEAYlIQhB5lEEgaDgAAFxFqf251RtRStHW2tb+9jS0nPseXpsa/scj63aY20dqlXRWkXAeSBa68Qg8xgZkwAJQ0LClOn6/bE3NFKGnZC9187e3/frxcusYa913QTzzVr3Wvdt7o6IiEikWgRdgIiINC8KDhERaRAFh4iINIiCQ0REGkTBISIiDaLgEBGRBlFwiETIzL5mZm/VW3Yz6xtkTUdjZnlmVniU7U+a2X/FsiZJDAoOkQi5+zPufn7QdYgETcEhIiINouCQhGBmPzKzIjOrMLNVZnZueH0LM5tqZl+Y2XYze8HMTgxv6xm+3TTZzDaZ2U4zu83MRpnZYjMrM7OH6p3jJjP78Ajn72BmT5lZqZltMLOfmlmL8LYUM/tvM9tmZuvM7I7weVua2f9nZvMPOdZdZvbKEc7T1cxmmtkOMysws1vrbWsdvv2008yWA6MO+exwM1sQ/jt6HmhVb1ummc0Ot3mHmf39QP0ih9I/DGn2zGwAcAcwyt3bARcA68ObvwNcAYwDugI7gYcPOcQZQD/gGuABYBowETgFuNrMxkVQxoNAB6B3+FxfByaHt90KXAScBowI13PATKCXmQ2qt+5G4KkjnGc6UBhuy1eBX5jZhPC2e4A+4T8XAN848CEzSwNmAE8DJwJ/Bb5S77g/CB83C8gGfgJoPCI5LAWHJIJaIB0YbGap7r7e3b8Ib7sNmObuhe6+H/gZ8FUza1nv8z93933u/hawG3jO3UvcvQj4OzD8aCc3sxTgWuDH7l7h7uuB/yYUAABXA78N17AT+NWBz4Zreh64IXysU4CewOzDnKc7cBbwo3C9C4E/EQqpA+e51913uPsm4Hf1Pj4aSAUecPdqd38RmFtvezVwEtAjvP3vroHs5AgUHNLsuXsB8D1CoVBiZtPNrGt4cw/g5fAtmDJgBaGgya53iK31vt57mOWMY5SQSeiH8oZ66zYA3cJfdwU21dtW/2uAPwPXm5kRCpsXwoFyqK7ADneviPA89evpChQdEgb1t/8GKADeMrO1Zjb1MOcXARQckiDc/Vl3P5tQUDhwX3jTJuAid+9Y70+r8NVEU9lG6Df2HvXWnQwcOMdmIKfetu6H1P4JUAWcA1xP6HbS4RQDJ5pZu6Ocp/sh26i3rVs4nP5le/hK6Qfu3hu4DLjrQD+RyKEUHNLsmdkAM5tgZunAPkJXCXXhzY8A95pZj/C+WWZ2eVOe391rgRfC52kXPtddwF/Cu7wAfNfMuplZR+BHhznMU8BDQLW7H7YDPnz76SPgl2bWysyGArcccp4fm9kJZpZDqH/ngI+BGuBOM0s1s6uA0w9sNLNLzaxvOFjKCV2V1SFyGAoOSQTphPoNtgFbgM7Aj8PbfkuoA/otM6sAPiHUGd7UvkOof2Qt8CHwLPB4eNsfgbeAxcDnwGuEfojX1vv808AQ/hkCR3IdoT6QYuBl4B53fye87T8J3X5aFz7fwSsXd68CrgJuAnYQehDgpXrH7Qe8A1QSCpnfu/ucCNotScjU/yUSW2Z2EfCIu/eot641UAKMcPc1gRUnEgFdcYhEWfj9iovD7210I/TY7MuH7HY7MFehIc2BrjhEoszM2gDvAwMJ9b+8CnzX3XeFt68HDLjC3T8Pqk6RSCk4RESkQXSrSkREGqTlsXdpfjIzM71nz56N+uzu3btp27Zt0xYU55KxzZCc7U7GNkNytruhbZ4/f/42d8+KZN+EDI6ePXsyb968Rn02Pz+fvLy8pi0oziVjmyE5252MbYbkbHdD22xmG469V4huVYmISIMoOEREpEEUHCIi0iAJ2cdxONXV1RQWFrJv376j7tehQwdWrFgRo6qio1WrVuTk5JCamhp0KSKSgJImOAoLC2nXrh09e/bkywOEfllFRQXt2rU74vZ45+5s376dwsJCevXqFXQ5IpKA4v5WlZldYWZ/NLPnzez8xh5n3759dOrU6aihkQjMjE6dOh3zykpEpLECCQ4ze9zMSsxs6SHrLwzPF11wYCIZd5/h7rcSmsntmuM87/F8vNlIlnaKSDCCuuJ4Eriw/orw9JsPE5qbeTBwnZkNrrfLT/nXuaJFRASYu34HD88pYPf+mqifK7CxqsysJzDb3YeEl8cAP3P3C8LLB+ZT+FX4z9v15h043PGmAFMAsrOzR06fPv1L2zt06EDfvn2PWVdtbS0pKSkNbU6T+MUvfkFGRgZ33nnnl9Zv2LCBq6++mk8//TTiYxUUFFBeXh7RvpWVlWRkHGt21MSTjO1OxjZDcrT7uZX7eXdjDY9MbEPLFtbgNo8fP36+u+dGsm88dY5348vzJRcSmnDnO8BEoIOZ9XX3Rw73YXd/FHgUIDc31w99Y3LFihURdXoH2Tmenp5Oenr6v5w/IyODFi1aNKiuVq1aMXz48Ij2Tca3aiE5252MbYbkaPfvV33MkG51TJxwFhDdNsd957i7/87dR7r7bUcKjQPMbJKZPRrpb9qxdv/99zNkyBCGDBnCAw88AMC9995L//79Ofvss1m1atXBfefPn8+wYcMYNmwYDz/8zzt0y5Yt4/TTT+e0005j6NChrFmj6RtEkl1tnbOsqJyhOR1icr54uuIoArrXW84Jr4uYu88CZuXm5t56tP3+c9YylhfvOuy2xt6qGty1PfdMOuWI2+fPn88TTzzBp59+irtzxhlncM455zB9+nQWLlxITU0NI0aMYOTIkQBMnjyZhx56iLFjx3L33XcfPM4jjzzCd7/7Xb72ta9RVVVFbW3tkU4pIklibWklu6tqGZrTMSbni6crjrlAPzPrZWZpwLWE5opOCB9++CFXXnklbdu2JSMjg6uuuopXX32VK6+8kjZt2tC+fXsuu+wyAMrKyigrK2Ps2LEA3HjjjQePM2bMGH7xi19w3333sWHDBlq3bh1Ie0QkfiwuDN1lSegrDjN7DsgDMs2sELjH3R8zszuAN4EU4HF3X9bA404CJh2rE/xoVwbx/gLg9ddfzxlnnMGrr77KxRdfzB/+8AcmTJgQdFkiEqDFhWW0SUuhT1ZsHgAI5IrD3a9z95PcPdXdc9z9sfD619y9v7v3cfd7G3HcWe4+pUOH2KRuQ5xzzjnMmDGDPXv2sHv3bl5++WUuueQSZsyYwd69e6moqGDWrFkAdOzYkY4dO/Lhhx8C8Mwzzxw8ztq1a+nduzd33nknl19+OYsXLw6kPSISPxYXlTOkawdSWsTmHa546uM4bpFecQRhxIgR3HTTTZx++ukAfPOb32TkyJFcc801DBs2jM6dOzNq1KiD+z/xxBPcfPPNmBnnn//PF+ZfeOEFnn76aVJTU+nSpQs/+clPYt4WEYkf1bV1LC/exY2je8TsnAkVHJF2jgflrrvu4q677vrSumnTpjFt2rR/2XfkyJEsWrTo4PKvf/1rAKZOncrUqVOjW6iINBurt1awv6aOU2PUvwHx1TkuIiINtORgx3hsnqgCBYeISLO2qLCcdq1a0rNTm5idM6GC41gvAAY1vEqsJUs7RQSWFJUxNKdDTAc3TajgONpTVa1atWL79u0J/0P1wHwcrVq1CroUEYmyfdW1rNxcEdPbVJBgneNHk5OTQ2FhIaWlpUfdb9++fc3+h+6BGQBFJLGt3FJBTZ0ztFtsX0FImuBITU2NaEa8/Pz8iAcHFBEJ0uLCMgCGdo/tFUdC3aqK90EORUSa0idrt5PVLp2uHWJ7lyShgiOe3xwXEWlKlftreG9lCRee0iXms34mVHCIiCSLd5ZvZV91HZed1jXm51ZwiIg0QzMXFdO1QytGnnxCzM+t4BARaWZ27q7ig9WlXDqsKy1iNLBhfQkVHOocF5Fk8MayLdTUOZcNi/1tKkiw4FDnuIgkg5kLi+md2ZZTurYP5PwJFRwiIolu6659fLJuO5cO6xrzp6kOUHCIiDQjry7ejDuB3aYCBYeISLNRV+c899lGTunanr6dYzNN7OEoOEREmom3lm9hTUklU8b2DrSOhAoOPVUlIonK3XnwvQJ6dmrDpUODu00FCRYceqpKRBJV/upSlhXv4lt5fUkJ4N2N+hIqOEREEpG78+C7a+jWsTVXDO8WdDkKDhGRePfx2u0s2FjGv43rTVrL4H9sB1+BiIgckbvzwDtryGqXztW53YMuB1BwiIjEtRkLi/hs3Q6+N7EfrVJTgi4HUHCIiMSt8j3V/NfsFQw/uSPXjTo56HIOSpqpY0VEmpv73lxJ2d5qnr7i1EBGwT0SXXGIiMShBRt38txnG7npzJ4MDmgwwyNJqODQC4Aikgj2Vdfyk5eWkN2uFd8/r3/Q5fyLhAoOvQAoIongv15dzsotFfziqiFkpMdfj0JCBYeISHM3e3Exf/lkI7ee04sJA7ODLuewFBwiInFiw/bdTP3bEoaf3JEfXjgw6HKOSMEhIhIH9lbV8u1nF5DSwnjwuuGkpsTvj+f4u3kmIpJkauucO6d/zrLiXfzp67nknNAm6JKOKn4jTUQkCbg7P5+9nLeXb+WeSwdz7qD47NeoT8EhIhKgxz5cx5MfreeWs3tx01m9gi4nIgoOEZGAvDB3E/e+toKLhnRh2sWDgi4nYgoOEZEAvDBvEz96aTFn983kf645La6GFDkWBYeISIy9MG8TP/pbKDT++PXcuBn1NlJxHxxm1tvMHjOzF4OuRUTkeP35o/XNOjQgoOAws8fNrMTMlh6y/kIzW2VmBWY2FcDd17r7LUHUKSLSVNyd+95YyT0zl3HuwOxmGxoQ3BXHk8CF9VeYWQrwMHARMBi4zswGx740EZGmVV1bxw/+uoj/zf+C604/mUduGNFsQwPA3D2YE5v1BGa7+5Dw8hjgZ+5+QXj5xwDu/svw8ovu/tWjHG8KMAUgOzt75PTp0xtVV2VlJRkZGY36bHOVjG2G5Gx3MrYZgm33rv3Owwv3sWpnHVf2TeWyPqmYRb8jvKFtHj9+/Hx3z41k33h6c7wbsKneciFwhpl1Au4FhpvZjw8EyaHc/VHgUYDc3FzPy8trVBH5+fk09rPNVTK2GZKz3cnYZgiu3UuLyvnJU/PYvhseuOY0rhjeLWbnjmab4yk4DsvdtwO3BV2HiEhDvDi/kGkvL6FT2zRevO1MTs1JnOke4ik4ioDu9ZZzwusiZmaTgEl9+/ZtyrpERCK2e38N/z5jKS99XsTo3ify0PUjyMxID7qsJhVPj+POBfqZWS8zSwOuBWY25ACayElEgrS0qJxJD37IjIVFfG9iP5755uiECw0I6IrDzJ4D8oBMMysE7nH3x8zsDuBNIAV43N2XNfC4uuIQkZirqa3jf/O/4LfvruHEtmk8883RjOnTKeiyoiaQ4HD3646w/jXgteM47ixgVm5u7q2NPYaISEMUlFTyg78uYtGmMiYN68r/vewUTmibFnRZURVPfRwiIs1GVU3oKuPhOQW0SU/hoeuHc+nQrkGXFRMJFRy6VSUisTBv/Q5+/NIS1pRUMmlYV/7j0sFktUu8vowjSajg0K0qEYmmkop9/Or1lby0oIhuHVvzxE2jGD+wc9BlxVxCBYeISDRU1dTx1Mfr+e07a9hXU8u38vpwx4S+tElLzh+hCdVq3aoSkabk7ry5bCu/en0F67fvYVz/LO6ZNJjeWck3bEt9CRUculUlIk1l/oYd3Pf6Kj5bv4N+nTN4YvIo8vpnxWScqXiXUMEhInK8Vm2p4DdvruKdFVvJzEjnv64YwrWjutMyJZ7elw6WgkNEBCgoqeC37xYwe3ExGWktufuCAUw+q2fS9mMcTUL9jaiPQ0Qaas3WCh6eU8Ari4ppnZrCbeP6MOWc3gn/Et/xSKjgUB+HiERqaVE5D88p4I1lW2jVMoUpY3sz5ZzedErAsaWaWkIFh4jI0bg7H3+xnUc+WMsHq0tpl96SO8b3ZfJZvThRVxgRU3CISMKrqXNeWVjEn/6+jiVF5WRmpHP3BQO4cUwP2rdKDbq8ZkfBISIJq3xPNdPnbuQPH+xlx76F9M5syy+vOpUrh3dr1nN+By2hgkOd4yICoUdqn/xoPTM+L2JvdS2DTmzBb64ZyfgBnWnRQu9hHK+ECg51joskr6qaOt5avoWnP97Ap+t2kN6yBVec1o1vnNmTktULyBuUHXSJCSOhgkNEks+mHXt4fu4mps/dxLbK/XQ/sTVTLxrINbndDz5SW7I64CITjIJDRJqdqpo63lu5lec+28QHa0oxIG9AZ24c3YOx/bNI0e2oqFJwiEizUVBSwfNzN/HSgiK2766iS/tW3DmhH1eP6k63jq2DLi9pKDhEJK6V761m1qJiXpxfyMJNZbRsYZw7qDPXjOrO2H5ZGkMqAAkVHHqqSiQxVNfW8fc1pfxtQRFvL99KVU0dA7Lb8dNLBnH5ad2Sara9eJRQwaGnqkSaL3dnUWE5Mz4vYvbiYrZVVnFi2zSuG9Wdr47szpBu7TWkeZxIqOAQkeanoKSSmYuKmbWomHXbdpPWsgXnDuzMVSNyGNc/i7SWuhUVbxQcIhJzhTv3MHvxZmYtKmZZ8S7MYHSvTtw+rg8XntpFw4DEOQWHiMTE5vK9vLZkC68uLmbBxjIAhuV04N8vHcylQ08iu32rgCuUSCk4RCRqisv28sbSLby2ZDPzNuwEYNBJ7fnhhQO49NSunNypTcAVSmMoOESkSW3asYfXl27m9aVb+Dx8ZTGwSzv+z/n9ufjUk+idlRFwhXK8FBwiclzcnYKSSt5ctoXXl25hWfEuAE7p2p67LxjARUO6KCwSjIJDRBqsrs5ZVFjGm8u28tayLazdthuAESd3ZNrFg7hwSBe6n6jbUIkqoYJDLwCKRM/+mlo+/mI7by/fytvLt1JSsZ+WLYzRvTsx+exenD84Wx3cSSKhgkMvAIo0rbI9VcxZVcI7y0t4f3UplftraJOWwth+WVwwJJsJA7Lp0EaPziabhAoOETl+a0sreXdFCe+s2Mq8DTuprXOy2qUzaVhXzhvcmTP7ZGr2vCSn4BBJcjW1dczbsJP3VobCYm1pqL9iYJd23D6uDxMHZzO0WwfNnCcHKThEktCO3VW8v7qE91aW8v6qEnbtqyE1JdRf8fXRPTh3ULY6t+WIFBwiScDdWb55F7O+qOLBFR/x+cad1DlkZqRx/ildmDioM2f3yyIjXT8S5Ngi+ldiZn2AQnffb2Z5wFDgKXcvi2ZxItJ4lftr+EfBNuasLGHOqhK27toPwNCcOr4zoR8TBnbmVN2CkkaI9NeLvwG5ZtYXeBR4BXgWuDhahYlIw7g7X5RWkr+qlDmrSvhs3Q6qa5126S05p38meQM6k7Z9DVdccHbQpUozF2lw1Ll7jZldCTzo7g+a2efRLExEjm33/ho++mI7+atKyF9VSlHZXgD6Z2dw81m9GDcgi1E9TyQ1PEtefv4XQZYrCSLS4Kg2s+uAbwCTwuv08LZIjLk7a0oqeX9VKfmrS5i7bidVtXW0TUvhzL6ZfGt8H8b1zyLnBHVsS/REGhyTgduAe919nZn1Ap6OXlkickDFvurwVUUpH6z+8lXFTWf1JK9/FiN7nkB6S71bIbERUXC4+3LgznrL64D7olWUSDI78ATU+6tLeX9VKfM37KSmzslIb8lZfTtxx4S+jO2fRbeOrYMuVZLUUYPDzJYAfqTt7j60ySv61xraAr8HqoB8d38m2ucUibWdu6v4e8E23l9VygdrSimtCD0BNeik9tw6tjfj+mcxsscJB/sqRIJ0rCuOS8P//Xb4vwduT93AUQLlWMzs8fCxS9x9SL31FwK/BVKAP7n7r4CrgBfdfZaZPQ8oOKTZq61zFm4q4/3VodtPiwrLcIeObVI5u2/oCaix/TLprEEDJQ4dNTjcfQOAmZ3n7sPrbfqRmS0ApjbyvE8CDwFPHVhhZinAw8B5QCEw18xmAjnAkvButY08n0jgtu7aF7r9tLqUD9dso3xvNS0MhnXvyJ0T+pE3IIuhOR1J0XsVEufM/dgXDma2EPi2u/8jvHwm8Ht3P63RJzbrCcw+cMVhZmOAn7n7BeHlH4d3LQR2uvtsM5vu7tce4XhTgCkA2dnZI6dPn96ouiorK8nISK5JZ5KxzRD9dlfXOWt21rFkWy1LSmsorAz9v9Yx3Tg1M4UhmSmc0imFjLTYBYW+18mjoW0eP378fHfPjWTfSJ+qugV43Mw6AAbsBG6OuKLIdAM21VsuBM4Afgc8ZGaXALOO9GF3f5TQy4nk5uZ6Xl5eo4rIz8+nsZ9trpKxzRCddm/cvof3V4eGIP/oi+3sqaolNcUY1fNEbjgni7H9shh0UjvMgrmq0Pc6eUSzzZE+VTUfGBYODty9PCrVHP7cuwk9DnxMmshJYm1vVS2frN1+8BbUuvBMeCef2IavjMhhXP8sxvTpRFuNASUJ5FhPVd11hPUAuPv9TVhLEdC93nJOeF3ENJGTRFv9YT3eX13Kp+t2UFVTR6vUFozp3YlvjOnBuAGd6ZXZNuhSRaLmWL8GtQv/dwAwCpgZXp4EfNbEtcwF+oVfLiwCrgWub+JziDTYkYb16Ns5g6+P7nFwWA9NbiTJ4lhPVf0ngJl9AIxw94rw8s+AVxt7UjN7DsgDMs2sELjH3R8zszuANwk9jvu4uy9r4HF1q0qOm7tTUBK6qjh0WI8xfTK5Pa8PeQM0rIckr0hvvGYTegHvgKrwukZx9+uOsP414LXjOK5uVUmj7Kmq4aOC7eSvLmHOysMP65Hb80TSWuoFPJFIg+Mp4DMzezm8fCXw5+iUJBIbJXvqeOIf65izqpRP1m6nqqaONmkpnBUeLDBvQGcN6yFyGJE+VXWvmb0OXELojfGb3D3uhlXXrSo5mqqaOuau38F7K0uYs7KEtdv2AsvpndWWG0f3YPyAzozqpcECRY4l0hkA7wRuBV4i9B7Hn83sj+7+YDSLayjdqpJDbavcz5yVJby3soS/r9lG5f4a0lq2YHTvTozJqmbKpWfSo5OegBJpiEhvVX0TGB1+pwIzuw/4GIir4BBxd1ZuqeDdFVt5Z0XJwTGgstunM2lYVyYM7MxZfTvRJq0l+fn5Cg2RRog0OIwvjxNVG14XV3SrKjlV1dTx6brtvLM8FBYHOraH5XTg+xP7M2FgZ07p2j6wt7VFEk2kwfEE8Gm9zvErgMeiU1Lj6VZV8ijfU82cVSW8vWIr768qpXJ/Da1SW3BOvyzuPLcv4wd2pnM7jSwrEg2Rdo7fb2b5wIFZ7ifHY+e4JLbisr28tWwLby3fymfrdlBT52RmpDNp2ElMHJTNWX0z9RKeSAxEPICOuy8AFkSxFpF/UVBSwRtLt/Dmsq0sKQoNkdYnqy23ju3NeYOzOS2nIy00DLlITCXUyGvq42j+3J2lRbt4felm3li2hbWloUEDh5/ckakXDeS8wdn0yUqu4bFF4k1CBYf6OJqnujrn801lvL5kM68v3UJR2V5SWhhjendi8lm9OH9wNtmaCU8kbiRUcEjzEQqLncxevJk3lm5hc/k+0lJacHa/TL47sR/nDcrmhLZpQZcpIoeh4JCYcQ/Nsz178WZeW7I5FBYtWzCufxY/unAgEwZ1pn2r1KDLFJFjUHBIVLk7KzZXMGtxMbMWFVO4cy9pKS0Y2z+LqRcNZMLAzrRTWIg0KwkVHOocjx+bduzhlYVFvLKwmDUllaS0MM7qm8n3JvbnvMHZdGitsBBprhIqONQ5HqyyPVXMXryZGZ8XMW/DTgBG9TyBn18xhIuHdKFTRnrAFYpIU0io4JDYq66tY87KEl5aUMR7K0uoqq2jf3YGP7xwAJcN66rJjkQSkIJDGmV58S5enF/IKwuL2L67isyMNG4c04Mrh3fTuFAiCU7BIREr31vNzEXFvDB3E0uKyklNMc4bnM1XRuQwrn8WLVM0O55IMlBwyFG5O/M37OTZzzby6uLN7K+pY9BJ7bln0mCuOK2b3rUQSUIJFRx6qqrp7NpXzUvzC3n2s42s3lpJRnpLvjIyh+tGncyQbroVJZLMEio49FTV8VtevIunP9nAjM+L2Ftdy9CcDtz3lVO5dGhX2qYn1D8XEWkk/SQQauucVxdv5s8freez9TtoldqCy4Z15YbRPRia0zHo8kQkzig4klj53mqen7uRRz7Yy459C+h+YmumXTyIq3O706GNXtATkcNTcCShTTv28NiH63hh3ib2VNUy6MQW3Hf1SCYM7EyK5rYQkWNQcCSR5cW7+MMHXzB78WYMuGxYV245pxelqz8nb3B20OWJSDOh4EgCCzbu5OH3Cnh3ZQlt01K4+ayeTD6rF107tgYgf3XABYpIs6LgSGBz1+/ggXdW84+C7ZzQJpUfnNefr4/pqf4LETkuCo4ENH/DTv7n7dV8WLCNzIx0pl08iOvPOFmP04pIk0ionyTJ/gLgis27+H9vruLdlSV0apvGTy8ZxNfO6EHrtJSgSxORBJJQwZGsLwAW7tzDf7+1mhkLi2iX3pK7LxjA5LN60iYtob69IhIn9JOlGavYV83v87/gsQ/XYcC/je3D7eP6qA9DRKJKwdEM1dU5f52/iV+/sYrtu6u4ang3/s8FAw4+JSUiEk0KjmZmcWEZ//7KMhZtKiO3xwk8MXmUhgURkZhScDQTu/ZV85s3VvGXTzfQqW069189jCuHd9MotSIScwqOZuDNZVv4j1eWUlKxn2+M6ckPzu9Pu1bqxxCRYCg44tiO3VX8+4ylvLpkMwO7tOPRG3MZ1l23pUQkWAqOOPXO8q1MfWkJ5XuruPuCAUwZ25tUTc0qInFAwRFn9lbV8p+zljF97iYGdmnH07eczqCT2gddlojIQQqOOLJ6awXffmYBBaWV3J7Xh+9P7E9aS11liEh8ifvgMLPewDSgg7t/Neh6ouWFeZv4j1eWkpHekqduPp1z+mUFXZKIyGFF9ddZM3vczErMbOkh6y80s1VmVmBmU492DHdf6+4G6v+EAAALrElEQVS3RLPOIFXX1vEfryzlhy8uZmSPE3jtu+coNEQkrkX7iuNJ4CHgqQMrzCwFeBg4DygE5prZTCAF+OUhn7/Z3UuiXGNgtlfu51vPLODTdTuYMrY3P7xgAC3VAS4icS6qweHuH5hZz0NWnw4UuPtaADObDlzu7r8ELo1mPfHki9JKvvH4Z5RW7OeBa07jiuHdgi5JRCQi5u7RPUEoOGa7+5Dw8leBC939m+HlG4Ez3P2OI3y+E3AvoSuUP4UD5nD7TQGmAGRnZ4+cPn16o+qtrKwkIyOjUZ+NVMHOWv5nwT5SDL43shW9OwQ77Hks2hyPkrHdydhmSM52N7TN48ePn+/uuZHsG/ed4+6+Hbgtgv0eBR4FyM3N9by8vEadLz8/n8Z+NhLvLN/K/3t3Adnt2/DUzafTo1PbqJ0rUtFuc7xKxnYnY5shOdsdzTYHcUO9COhebzknvO64mdkkM3u0vLy8KQ7X5F5dvJl/+8t8+me342+3nxkXoSEi0lBBBMdcoJ+Z9TKzNOBaYGZTHNjdZ7n7lA4dOjTF4ZrUq4s3c+f0zxnevSPP3jqazIz0oEsSEWmUaD+O+xzwMTDAzArN7BZ3rwHuAN4EVgAvuPuyaNYRtNeW/DM0nrz5dDI097eINGPRfqrquiOsfw14ranPF49zjr+/upTvPPc5pyk0RCRBJNRLA/F2q2pZcTnf+st8+nXO4MnJoxQaIpIQEio44qlzvLhsLzc/OZf2rVN5cvLpmj9DRBJGQgVHvFxxVO6v4eYn57J7fy2P3zSKLh1aBVqPiEhTSqjgiAfuzk9eWsLqrRX87w0jNCS6iCQcBUcTmz53EzMXFfP9if01WKGIJKSECo6g+zhWbN7Fz2Yu45x+mXxrfPw82SUi0pQSKjiC7OPYvb+Gbz+7gPatU7n/6tNIaWExr0FEJBb0fGgTuf/t1azbtptnvnkGWe30VriIJK6EuuIIyrLicp78aD3XnX4yZ/bJDLocEZGoSqjgCKKPo67OmfbyUjq2TuVHFwyM2XlFRIKSUMERRB/Hc3M3snBTGT+9dBAd2uglPxFJfAkVHLFWWrGf+15fyZjenbjiNM3gJyLJQcFxHB6eU8Ceqlp+fsUQzPQUlYgkh4QKjlj2cZRW7Oe5zzZy5fBu9O2cXFNSikhyS6jgiGUfx58+XEt1bR235/WJ+rlEROJJQgVHrOzcXcVfPt7ApUO70jtLVxsiklwUHI3wxEfr2V1Vy7c1rIiIJCEFRwNV7KvmyX+s4/zB2Qzo0i7ockREYk7B0UDPfbaRXftquGOCrjZEJDkpOBropQVFDD+5I0NzOgZdiohIIBIqOKL9OO6arRWs3FLBZcO6RuX4IiLNQUIFR7Qfx521qJgWBpcMPSkqxxcRaQ4SKjiiyd2ZuaiYMX060bmd5hAXkeSl4IjQ0qJdrN++h0lDdZtKRJKbgiNCMxcVkZpiXDREt6lEJLkpOCJQV+fMXryZcf2zNHS6iCQ9BUcE5m3YyebyfUzS01QiIgqOSLy2ZDOtUlswcVB20KWIiAQuoYIjWu9xzNuwg5E9TqBtessmPa6ISHOUUMERjfc49lXXsmpLBad205viIiKQYMERDSu3VFBd6wzLid085iIi8UzBcQxLCssAOFXBISICKDiOaVFhOZ3aptGtY+ugSxERiQsKjmNYUljOqTkdMLOgSxERiQsKjqPYU1XDmpIKDaEuIlKPguMolhXvos5haDf1b4iIHKDgOIrFhaH3QYaqY1xE5CAFx1EsLiyjS/tWdG6vYdRFRA5QcBzFksJyXW2IiBxCwXEE5XurWbttt4JDROQQcT/4kpldAVwCtAcec/e3YnHeZUUH+jf0RJWISH1RveIws8fNrMTMlh6y/kIzW2VmBWY29WjHcPcZ7n4rcBtwTTTrrW9xODhO1RNVIiJfEu0rjieBh4CnDqwwsxTgYeA8oBCYa2YzgRTgl4d8/mZ3Lwl//dPw52JiSWE53U9szQlt02J1ShGRZsHcPbonMOsJzHb3IeHlMcDP3P2C8PKPAdz90NA48HkDfgW87e7vHOU8U4ApANnZ2SOnT5/eqHorKyvJyMjg13P3Ul0L00Yn/lAjB9qcbJKx3cnYZkjOdje0zePHj5/v7rmR7BtEH0c3YFO95ULgjKPs/x1gItDBzPq6+yOH28ndHwUeBcjNzfW8vLxGFZefn09eXh5/LPiE/dV15OWd2ajjNCcH2pxskrHdydhmSM52R7PNcd857u6/A34XdB0iIhISxOO4RUD3ess54XXHLVozAIqIyD8FERxzgX5m1svM0oBrgZlNceBozAAoIiJfFu3HcZ8DPgYGmFmhmd3i7jXAHcCbwArgBXdf1kTn0xWHiEiURbWPw92vO8L614DXonC+WcCs3NzcW5v62CIiEqIhR0REpEESKjh0q0pEJPoSKjjUOS4iEn1Rf3M8CGZWCmxo5MczgW1NWE5zkIxthuRsdzK2GZKz3Q1tcw93z4pkx4QMjuNhZvMife0+USRjmyE5252MbYbkbHc025xQt6pERCT6FBwiItIgCo5/9WjQBQQgGdsMydnuZGwzJGe7o9Zm9XGIiEiD6IpDREQaRMEhIiINkrTBcax5z80s3cyeD2//NDyTYbMWQZvvMrPlZrbYzN41sx5B1NnUIp3j3sy+YmZuZs3+sc1I2mxmV4e/38vM7NlY1xgNEfwbP9nM5pjZ5+F/5xcHUWdTMrPHzazEzJYeYbuZ2e/CfyeLzWzEcZ/U3ZPuD6H5zb8AegNpwCJg8CH7fAt4JPz1tcDzQdcdgzaPB9qEv769ubc50naH92sHfAB8AuQGXXcMvtf9gM+BE8LLnYOuO0btfhS4Pfz1YGB90HU3QbvHAiOApUfYfjHwOmDAaODT4z1nsl5xnA4UuPtad68CpgOXH7LP5cCfw1+/CJwbnv+8uTpmm919jrvvCS9+QmiSreYuku81wM+B+4B9sSwuSiJp863Aw+6+E8DdS2JcYzRE0m4H2oe/7gAUx7C+qHD3D4AdR9nlcuApD/kE6GhmJx3POZM1OA4373m3I+3joTlEyoFOMakuOiJpc323EPotpbk7ZrvDl+7d3f3VWBYWRZF8r/sD/c3sH2b2iZldGLPqoieSdv8MuMHMCglN7fCd2JQWqIb+v39McT/nuMSemd0A5ALjgq4l2sysBXA/cFPApcRaS0K3q/IIXVl+YGanuntZoFVF33XAk+7+32Y2BnjazIa4e13QhTUnyXrFEcm85wf3MbOWhC5rt8ekuuiIaK53M5sITAMuc/f9Maotmo7V7nbAECDfzNYTugc8s5l3kEfyvS4EZrp7tbuvA1YTCpLmLJJ23wK8AODuHwOtCA0GmMgi+n+/IZI1OCKZ93wm8I3w118F3vNwT1Mzdcw2m9lw4A+EQiMR7nnDMdrt7uXununuPd29J6G+ncvcfV4w5TaJSP59zyB0tYGZZRK6dbU2lkVGQSTt3gicC2BmgwgFR2lMq4y9mcDXw09XjQbK3X3z8RwwKW9VuXuNmR2Y9zwFeNzdl5nZ/wXmuftM4DFCl7EFhDqerg2u4uMXYZt/A2QAfw0/B7DR3S8LrOgmEGG7E0qEbX4TON/MlgO1wN3u3pyvqCNt9w+AP5rZ9wl1lN/UzH8hxMyeI/RLQGa47+YeIBXA3R8h1JdzMVAA7AEmH/c5m/nfmYiIxFiy3qoSEZFGUnCIiEiDKDhERKRBFBwiItIgCg4REWkQBYeIiDSIgkNERBpEwSESA2bW08xWmtkzZrbCzF40szZB1yXSGAoOkdgZAPze3QcBuwjN+SLS7Cg4RGJnk7v/I/z1X4CzgyxGpLEUHCKxc+j4PhrvR5olBYdI7JwcngMC4HrgwyCLEWksBYdI7KwCvm1mK4ATgP8NuB6RRknKYdVFAlLj7jcEXYTI8dIVh4iINIjm4xARkQbRFYeIiDSIgkNERBpEwSEiIg2i4BARkQZRcIiISIP8/4/h3SLvwY5JAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.semilogy(p, odds, label=\"odds\")\n",
    "plt.xlabel(\"p\")\n",
    "plt.ylabel(\"odds\")\n",
    "plt.legend()\n",
    "plt.title(\"semilogy odds\")\n",
    "plt.grid()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-02-25T01:47:42.412749Z",
     "start_time": "2019-02-25T01:47:42.222151Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEWCAYAAABxMXBSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xl81dWd//HXJ4Q97EsIBARkCfsWcKFiEFe0WLeqXd1KW8euU7uMbe105jdjbWemtdoqtWptHbFuIyqCC0RERVlk32QnIZAQDCSEkO3z++NebIwJuReS+73JfT8fjzy43+/35J7PSfR+cr7nfM8xd0dERCRSSUEHICIizYsSh4iIREWJQ0REoqLEISIiUVHiEBGRqChxiIhIVJQ4JCGZ2S4zuzDoOCJlZm5mQ+q5dpOZLY11TJK4lDhETpGZTTezxWZ22Mx21XF9YPh6qZltbk6JSuRklDhETt1R4BHgznquPwl8APQA7gKeMbNeMYpNpMkocUhCM7O2ZvZbM9sX/vqtmbWtcf2HZpYXvnZbzVtG7v6+u/8V2FHH+w4DJgJ3u/sxd38WWAdcE2G9d9ao95Za793DzOaZ2REzex84s8Y1M7P/MbP88PV1Zja6cX9qkuiUOCTR3QWcDYwHxgFTgJ8CmNmlwPeBC4EhQFYU7zsK2OHuxTXOrQmfj6TeHwAXAUPD9df0AFAGpAG3hL9OuBiYBgwDugCfBwqjiFukQUockui+CPzS3fPdvQD4V+DL4WufBx519w3uXgr8Ior3TQEO1zp3GOgURb3r3f1ozXrNrBWhXsvP3f2ou68H/lKjjopwHRmAufsmd8+LIm6RBilxSKLrC+yucbw7fO7Etb01rtV83ZASoHOtc52BEz2QaOqtWa4XkFzfdXdfBNxPqFeSb2ZzzKx2HCKnRYlDEt0+4IwaxwPC5wDygPQa1/pH8b4bgMFm1qnGuXHh85HU27/WtRMKgMqTXMfd73P3ScBIQres6hu8FzklShyS6J4EfmpmvcysJ/Bz4G/ha38HbjazEWbWAfhZzW80syQzawe0Dh1aOzNrA+DuW4HVwN3h81cBY4FnI6z3JjMbGa737hN1unsV8BzwCzPrYGYjga/WiGmymZ1lZq0JzfoqA6ob5SclEqbEIYnu34EVwFpCs55Whc/h7q8A9wGLgW3AsvD3HA//Ow04Bswn9Ff/MeDVGu99A5AJfATcA1wbHs+IpN7fAovC9S6qFfMdhMZQ9gOPAY/WuNYZ+FO4zt2EBsZ/HcXPQ6RBpo2cRCJjZiOA9UBbd68MOh6RoKjHIXISZnZV+JmLbsCvgBeVNCTRKXGInNzXgXxgO1AFfDPYcESCp1tVIiISFfU4REQkKslBB9AUevbs6QMHDoy4/NGjR+nYsWPTBRSn1O7EonYnlmjbvXLlyoPuHtEinC0ycQwcOJAVK1ZEXD47O5usrKymCyhOqd2JRe1OLNG228x2N1wqRLeqREQkKkocIiISFSUOERGJSosc46hLRUUFOTk5lJWVfepaly5d2LRpUwBRBasx2t2uXTvS09Np3bp1I0UlIvEuYRJHTk4OnTp1YuDAgZjZJ64VFxfTqVOner6z5Trddrs7hYWF5OTkMGjQoEaMTETiWWC3qsysv5ktNrONZrbBzL5TRxkzs/vMbJuZrTWziadaX1lZGT169PhU0pBTZ2b06NGjzl6ciLRcQfY4KoF/dvdV4T0LVprZa+6+sUaZywhtnTkUOAv4Y/jfU6Kk0fj0MxVJPIH1ONw9z91XhV8XA5uAfrWKXQk87iHLgK5mlhbjUEVE4t7rGw/w0Jvbqa5u+mWk4mKMw8wGAhOA92pd6scnt8jMCZ/71B7KZjYbmA2QmppKdnb2J6536dKF4uLi2t8GQFVVVb3XGlNaWhp5eae2/fNDDz3EH/7wB3bu3MnOnTvp0aMHEBpn+OEPf8irr75Khw4d+OMf/8j48eMjes+qqiq++MUvcumll/K5z33uE9feeust7rvvPp5++ukG36esrOxTP+94VlJS0qzibSxqd8v24Adl7DpSzXAPfWQ2ZbsDTxxmlkJoV7TvuvuRU30fd58DzAHIzMz02k9Mbtq0qd6B4FgOjp9qPTNmzODaa68lKyuLlJSUj99n/vz57N69m+3bt/Pee+/xne98h/feq51/61ZcXEzr1q1p3779p+Lq0KEDycnJEcXbrl07JkyYEH2jAqIniRNLorT7rmWLOGtIV7KyQkPBTdnuQJ/jCG9v+SzwhLs/V0eRXD65t3J6+Fyz5u7ceeedjB49mjFjxvDUU08BUF1dze23305GRgYXXXQRM2fO5JlnngFgwoQJ1LX+1gsvvMBXvvIVzIyzzz6boqIi8vLy6q3D3bnjjjsYPnw4s2bNIj8//+P3WrBgARkZGUycOJHnnvvHr+PNN99k/PjxjB8/ngkTJsSkdyYikSssOU5u0THGpneJSX2B9TgsNKr6Z2CTu/93PcXmAXeY2VxCg+KH3f3U7vXU8K8vbmDjvn90bqqqqmjVqtVpvefIvp25+7OjIir73HPPsXr1atasWcPBgweZPHky06ZN4+2332bXrl1s3LiR/Px8RowYwS233HLS98rNzaV//3/k1vT0dHJzc3nnnXfqrOPdd99ly5YtbNy4ke3btzNlyhRuueUWysrK+NrXvsaiRYsYMmQI119//cfv+Zvf/IYHHniAqVOnUlJSQrt27U7thyQiTWJtzmEAxqZ3jUl9QfY4pgJfBi4ws9Xhr5lm9g0z+0a4zHxgB6F9l/8E3B5QrI1q6dKl3HjjjbRq1YrU1FTOP/98li9fztKlS7nuuutISkqiT58+TJ8+vdHrWLJkycfn09LSuOCCCwDYvHkzgwYNYujQoZgZX/rSlz5+r6lTp/L973+f++67j6KiIpKTA7/DKSI1rMkpwgzGtPQeh7svBU46l9NDu0z9U2PXXbtn0JwfAOzXrx979/5j/kBOTg79+tWenHZ6fvzjH3P55Zczf/58pk6dysKFC8nIyGjUOkTk1K3NOcyZvVJIaRubj3StVRWA8847j6eeeoqqqioKCgpYsmQJU6ZMYerUqTz77LNUV1dz4MCBiGZEzJo1i8cffxx3Z9myZXTp0oW0tLR665g2bdrH5/fv38/ixYsByMjIYNeuXWzfvh2AJ5988uM6tm/fzpgxY/jRj37E5MmT2bx5c5P8XEQkeu7O2pyimI1vQBzMqkpEV111Fe+++y7jxo3DzLj33nvp06cP11xzDW+88QYjR46kf//+TJw4kS5dQv8x3Hfffdx7773s37+fsWPHMnPmTB5++GFmzpzJ/PnzGTJkCB06dODRRx89aR1XXXUVixYtYuTIkfTr149zzjkHCM2MmjNnDpdffjkdOnTgvPPO+3gQ/Le//S2LFy8mKSmJUaNGcdlllwXzgxORT9l3uIyDJeWMi9H4BhDKVi3ta9KkSV7bxo0bP3XuhCNHjtR7LdaKi4vd3f3gwYM+ePBgz8vLa7K6GqvdJ/vZxqPFixcHHUIg1O6Waf7afX7Gj17yVbsPfeJ8tO0GVniEn7HqccSZK664gqKiIsrLy/nZz35Gnz59gg5JROLYmpzDJCcZI9I6x6xOJY44kwhPuIpI41mbU0RGWifatT69RwqikVCD46HemDQm/UxFglNd7azLORyz5zdOSJjE0a5dOwoLC/VB14g8vB+HHggUCcbOwqMUH69kXAxnVEEC3apKT08nJyeHgoKCT10rKytLyA+/xmj3iR0ARST21uYUAbF7YvyEhEkcrVu3rneXuuzs7Ga1SF9jSdR2i7QUa/Yepl3rJIb2TolpvQlzq0pEpKVZtqOQsf26ktwqth/lShwiIs3QhweK2by/mMvGxH7KvhKHiEgzNG/NPpIMLh8b+01RlThERJoZd+eF1fs498ye9O4U+4k9ShwiIs3MmpzD7DlUyqzxfQOpX4lDRKSZeWF1Lm1aJXHJqGCWJAp669hHzCzfzNbXcz3LzA7X2Ojp57GOUUQknlRVOy+tzWN6Ri+6tG8dSAxBP8fxGHA/8PhJyrzl7lfEJhwRkfi2bEchBcXHuXJ8427YFo1AexzuvgQ4FGQMIiLNyQurc0lpm8wFGb0Di6E5jHGcY2ZrzOwVMxvVcHERkZapqLScl9fmccmoPjFdDbc2C3rRPzMbCLzk7qPruNYZqHb3EjObCfzO3YfW8z6zgdkAqampk+bOnRtxDCUlJaSkxPaR/XigdicWtbv5e+7DcuZtr+Dfp7YnvdPJ/+6Ptt3Tp09f6e6ZERWOdMenpvoCBgLrIyy7C+jZULm6dgA8mZa+Q1h91O7EonY3b4ePlfvouxf41x9fEVH5ptwBMK5vVZlZHzOz8OsphG6tFQYblYhI7P3l7V0Ul1VyxwVDgg4l2FlVZvYkkAX0NLMc4G6gNYC7PwhcC3zTzCqBY8AN4cwoIpIwSo5X8ue3dzIjozej+8V27426BJo43P3GBq7fT2i6rohIwvrbst0UlVbwrRl1DvHGXFzfqhIRSXRHyir405IdnDe0J+P7x3bDpvoocYiIxLFfL9jCR6Xl/PCSjKBD+ZgSh4hInFq9t4i/vbebr5wzkDEx3lf8ZJQ4RETiUGVVNT95bh29O7Xlny8eFnQ4nxD0WlUiIlKHR9/exaa8Izz4pYl0ahfMYob1UY9DRCTO7Dp4lP95fSszMnoHtnT6yShxiIjEkbKKKm5/YhVtkpP4t8+NJvwMdFzRrSoRkTjyby9tZGPeER65KZO+XdsHHU6d1OMQEYkTL6zO5Yn39vD18wdzQUZq0OHUS4lDRCQOfHigmH95bh2TzujGDy4eHnQ4J6XEISISsANHyrjp0eV0aJvM72+cQOtW8f3RHN/RiYi0cMVlFdz06HKKSst59KbJcTuuUZMGx0VEAlJRVc3tT6xi64FiHrlpclysfBsJ9ThERAJQUVXNd+eu5q0PD/KfV4/h/GG9gg4pYupxiIjEWEVVNd9+8gNeWb+fn14+gs9n9g86pKgocYiIxFB5ZTXfenIVCzcc4OdXjOSWzwwKOqSoBXqrysweMbN8M1tfz3Uzs/vMbJuZrTWzibGOUUSksRw9Xsnsv65g4YYD3P3Z5pk0IPgxjseAS09y/TJgaPhrNvDHGMQkItLoCoqPc8OcZSzZWsB/XDWGm6c2z6QBwW8du8TMBp6kyJXA4+F9xpeZWVczS3P3vJgEKCLSCHYUlPDVR9/nYHE5f/pKJjNGxO9T4ZGw0GdygAGEEsdL7j66jmsvAfe4+9Lw8RvAj9x9RR1lZxPqlZCamjpp7ty5EcdQUlJCSkrKKcXfnKndiUXtDsbagkr+uOY4yUnwvYntGNy1VUzqjbbd06dPX+numZGUbTGD4+4+B5gDkJmZ6VlZWRF/b3Z2NtGUbynU7sSidseWu/OH7O38z6otZPTpzJwvT6J/9w4xq78p2x3viSMXqDlPLT18TkQkbh0ureCHz65h4YYDzBrXl19dM5b2bWLT04iFeE8c84A7zGwucBZwWOMbIhLPVu4+xLefXM2BI2X89PIR3PqZQXG5p8bpCDRxmNmTQBbQ08xygLuB1gDu/iAwH5gJbANKgZuDiVRE5OSqqp0H39zOf7+2lb5d2/HMN89lfP+uQYfVJIKeVXVjA9cd+KcYhSMickq2F5Twg6fX8MGeIq4Ym8Z/XD2GznG2T3hjivdbVSIicauq2nnsnV3cu2Az7Vq34nc3jGfWuL4t7tZUbUocIiKnYOO+I/zk+XWs2VvEjIze/OfVY+jduV3QYcWEEoeISBRKyyv53Rsf8vBbO+navnXC9DJqUuIQEYmAu/Pi2jz+4+VN7D9SxvWZ/fnJzAy6dmgTdGgxp8QhItKADfsO88sXN/LezkOM6tuZ+78wgcyB3YMOKzBKHCIi9cg7fIzfLNzKcx/k0LV9a/7fVaO5YfIAWiUlzm2puihxiIjUcri0ggeXbOeRpTtxh9nTBnN71hC6tG+5U2yjocQhIhJWWl7Jo2/v4qE3t1N8vJJZ4/ryg4uHx3SNqeZAiUNEEl5peSVPLNvDQ0u2c7CknBkZvfnBJcMZkdY56NDikhKHiCSskuOVPLFsN396awcHS8o5b2hPvnvhUCadkbgD35FQ4hCRhPPR0XIefWcXf3lnF4ePVXDe0J58Z8bQhJ4pFQ0lDhFJGHsPlfLwWzv4+4ocjlVUccmoVG7PGsK4FroYYVNR4hCRFu+DPR/x56U7mb8uj1ZJxqxx/fj6+YMZltop6NCaJSUOEWmRKqqqWbB+P7999xjbF7xDp7bJfO28wdw8dRB9uiTGmlJNRYlDRFqU/OIy5r6/lyfe282BI8dJ7WD866xRXDMpnZS2+shrDPopikiz5+68v/MQf3tvDwvW51FR5Uwb1ov/97kzSDqwkQvOHRh0iC1K0DsAXgr8DmgFPOzu99S6fhPwa/6xz/j97v5wTIMUkbhVVFrOc6tyefL9PXyYX0Lndsl86ewz+PLZZzC4VwoA2fmbAo6y5QkscZhZK+AB4CIgB1huZvPcfWOtok+5+x0xD1BE4lJ1tfPezkM8tXwP89fvp7yymnH9u3LvtWP57Ni+tG/TKugQW7wgexxTgG3uvgPAzOYCVwK1E4eICPuKjvHsyhyeXpnDnkOldGqXzI2T+3PDlAF6wjvGLLStdwAVm10LXOrut4WPvwycVbN3Eb5V9Z9AAbAV+J67763n/WYDswFSU1MnzZ07N+JYSkpKSElJOcWWNF9qd2Jpju0+XuWsPFDF27kVbCysxoER3ZM4L701k1Jb0bZVw6vUNsd2N4Zo2z19+vSV7p4ZSdl4Hxx/EXjS3Y+b2deBvwAX1FXQ3ecAcwAyMzM9Kysr4kqys7OJpnxLoXYnlubS7qpqZ9mOQp5blcuC9XkcLa+if/f2fHtGOtdMTGdAj+gWHGwu7W5sTdnuiBKHmXUEjrl7tZkNAzKAV9y94jTqzgX61zhO5x+D4AC4e2GNw4eBe0+jPhGJU+7OxrwjvLB6H/NW72P/kTI6tU3mirF9uWpiP6YM7E5Sgu+BEU8i7XEsAc4zs27Aq8By4Hrgi6dR93JgqJkNIpQwbgC+ULOAmaW5e174cBag6REiLcjOg0d5cc0+Xlyzjw/zS0hOMrKG9+KnV4zgwhGptGutge54FGniMHcvNbNbgT+4+71mtvp0Knb3SjO7A1hIaDruI+6+wcx+Caxw93nAt81sFlAJHAJuOp06RSR4OR+V8vLaPF5am8e63MMATBnYnX//3GguH5NGt46Jt4d3cxNx4jCzcwj1MG4NnzvtPwXcfT4wv9a5n9d4/RPgJ6dbj4gEK7foGK+sy+PldXl8sKcIgLHpXfjp5SO4fGwaaV3aBxyhRCPSxPFdQh/gz4d7BYOBxU0Xlog0d3sPlbJg/X5eXpfH6r2hZDEyrTM/vHQ4V4zpG/Ugt8SPiBKHu78JvFnjeAfw7aYKSkSap235JSzcsJ9X1uexPvcIAKP7hZLFzNFpDOzZMeAIpTGcNHGY2YtAvQ96uPusRo9IRJoNd2d97hEWbtjPgg372ZZfAsD4/l35l5kZXDY6Tft1t0AN9Th+E/73aqAP8Lfw8Y3AgaYKSkTiV0VVNe/vPMSrG/bz6sYD5B0uI8lgyqDufPnsUVwyqo+WLW/hTpo4wreoMLP/qvVE4YtmtqJJIxORuHGkrIIlWwt4beMBFm/O50hZJW2Tk5g2rBffv2gYF45I1WyoBBLp4HhHMxtcY12pQYBuVoq0YHsPlfLGpgO8sTmfZTsKqahyunVozcWj+nDxyFTOG9pLCwomqEgTx/eAbDPbARhwBvD1JotKRGKuqtpZvbeIRZsP8MamfDbvLwZgcM+O3DJ1EBeOTGXigG600hPcCS/SWVULzGwooaVGADa7+/GmC0tEYuFIWQVvbT3Ios35ZG/Jp/BoOa2SjMwzunHXzBHMGNH7430tRE5oaFbV1fVcOtPMcPfnmiAmEWki7s72ghIWbc5n0eZ8Vuz6iMpqp0v71mQN78WMEamcP7QXXTq0DjpUiWMN9Tg+G/63N3Au8AahW1XTgXcAJQ6ROHesvIplOwr528bj/Oz9xew9dAyA4amduO28wcwY0ZsJ/buS3Cop4EiluWhoVtXNAGb2KjDyxIKDZpYGPNbk0YnIKdldeJTFm/NZvKWAZTsKOV5ZTZtWMG1Yd2ZPO5MLMnrTr6uW+ZBTE+ngeP8aq9RC6BmOAU0Qj4icgrKKUK8ie0sBb24tYOfBo0BoYPsLZw1g+vDelO1dz8UzJgccqbQEkSaON8xsIfBk+Ph64PWmCUlEGuLu7Dx49ONEcaJX0TY5iXPO7MFXzzmDrOG9P7HER/Y+zYaSxhHprKo7zOwqYFr41Bx3f77pwhKR2kqOV/LOtoMs+TCULE6MVQzuFepVnD+sF2cP7qE9LKTJRbN17DuE9sVw4P2mCUdETjixK96bWwtYsrWAlbs/oqLK6dCmFeee2ZPZ087k/KG9tMqsxFykW8d+Hvg1kE1oVtXvzexOd3+mCWMTSTiHjpbzVrhH8daHBykoDj0uldGnE7d+ZjDThvUk84zutEnWDCgJTqQ9jruAye6eD2BmvQiNcZxW4jCzS4HfEdoU6mF3v6fW9bbA48AkoBC43t13nU6dIvGkoqqaD/YUsWRrAUs+LGBd7mHcoWuH1nxmSE+yhvdm2tCe9O6sRQMlfkSaOJJOJI2wQuC0/uQxs1bAA8BFQA6w3MzmufvGGsVuBT5y9yFmdgPwK0ID8yLN1t5DpSz5MHT76Z1thRQfryTJYMKAbnx3xjDOH96LMf26aGkPiVuRJo4Fdcyqmn+S8pGYAmyrsXDiXOBKoGbiuBL4Rfj1M8D9ZmbuXu8eISLxprS8kmU7Clmy9SBLthawIzxVtl/X9lwxLo1pQ3tx7pCedGmvp7WlebBIP4PDy498Jnz41unOqjKza4FL3f228PGXgbPc/Y4aZdaHy+SEj7eHyxys4/1mA7MBUlNTJ82dOzfiWEpKSkhJSbz1eNTupuHu7C2uZv3BKtYdrOLDj6qpdGiTBBndWzG6Z+grraNhFrtehX7fiSXadk+fPn1lre0z6hXp4HgqsCv8levucbeJk7vPAeYAZGZmelZWVsTfm52dTTTlWwq1u/HUN6g9PLUTt5zXi2lDe5E5sFugU2X1+04sTdnuhhY5HA88CHQBcsOn082sCLjd3VedRt25QP8ax+k16qhdJsfMksNxFJ5GnSKNorKqmtV7iz6eKru21qD2tGGhZKGd8KQlaqjH8RjwdXd/r+ZJMzsbeBQYdxp1LweGhjeFygVuAL5Qq8w84KvAu8C1wCKNb0hQ9hUdY8nWUK9i6baDFJdpUFsSU0OJo2PtpAHg7svM7LR2AHT3SjO7A1hIaDruI+6+wcx+Caxw93nAn4G/mtk24BCh5CISE2UVVSzfdejjZLH1QAkAfTq347LRfcga3pupZ/bUEuSScBpKHK+Y2cuEnqXYGz7XH/gKsOB0K3f3+dSaneXuP6/xugy47nTrEYnU7sJ/rP/07vZCjlVU0aZVEpMHdeO6Sf05f3gvhvZOiemgtki8aWhZ9W+b2WWEpsX2C5/OBR4If+iLNGv1rSp7Ro8OXJeZTtbw0PpPHdpEszqPSMvW4P8N7v4K8EoMYhGJifzSav7yzi6yt+TzzvaGV5UVkU865T+jzGyOu89uzGBEmkJ5ZTXv7zzE4i35LN6Sz46CY8AGBvbowI1TBnzcq9CqsiKRaWg6bvf6LgEzGz8ckcaRf6SMxVtC+2ov/fAgR8uraJOcxNmDe3BOjwpuu+JcBqlXIXJKGupxFAC7CSWKEzx83LupghKJVnW1s37fYV7flM/izfmsyz0MQFqXdlw5oR8XDO/NuUNCYxXZ2dlKGiKnoaHEsQOY4e57al8ws711lBeJmdLySpZ+eJBFm/N5Y3M+BcXHP36u4s5LhnNBRm8y+nTSDCiRRtZQ4vgt0A34VOIA7m38cERO7sCRMt7YlM/rmw7w9raDHK+sJqVtMucP68WMEb3JGt6b7h3bBB2mSIvW0HTcB05y7feNH47IJ7k7Ww+U8NrG/by28QBrckK3oNK7tefGKQO4aGQqkwdqYyORWIp0kcOr6zh9GFhXa58OkdNWVe2s2vMRr27Yz6sbD7C7sBSAcf27cuclw7lwRCrDUvUQnkhQIp2OeytwDrA4fJwFrAQGmdkv3f2vTRCbJJDyymre2X6QhRtCPYuDJeW0bmWcc2ZPZk8bzIUjUknVLngicSHSxJEMjDixnHp4mfXHgbOAJYASh0StrKKK7C0FLFifxxub8ykuq6Rjm1ZMz+jNJaP6kDW8F53aaR0okXgTaeLoX2sPjvzwuUNmVtEEcUkLVVpeyeLNBcxfl8eizfkcq6iia4fWXDqqD5eN6cO5Z/bUg3gicS7SxJFtZi8BT4ePrw2f6wgUNUlk0mIcK69i0eZ8Xl63j0Wb8ymrqKZnShuuntiPy0ancdbg7rRupcFtkeYi0sTxT0DNrWP/Ajwb3htjelMEJs3b8coq3txSwItr83hj0wFKy6vomdKW6yb1Z+aYNKYM6q59K0SaqYgSh7u7mS0Fygk9Of6+NlSS2qqqnXe3F/LC6lwWbNhPcVkl3Tq05nMT+nHF2DTOGtRDyUKkBYh0Ou7ngV8D2YSWG/m9md3p7s80YWzSDLg763OP8PwHuby4dh8Fxcfp1DaZi0f14bPj0pg6pKduQ4m0MJHeqroLmHzimQ0z6wW8DpxS4ggvnvgUMBDYBXze3T+qo1wVsC58uMfdZ51KfdL49hUd4/kPcnn+g1y25ZfQplUS0zN6ceX4flyQ0VsD3CItWKSJI6nWg36FwOn8Gflj4A13v8fMfhw+/lEd5Y65+/jTqEcaUWl5JQvW7+fZVTm8s70Qd5g8sBv/cdUYLh+Tpi1URRJEpIljgZktBJ4MH19PrS1fo3QloYcIITTQnk3diUMC5u6s2lPE0yv28tLaPEqOV9K/e3u+M2MoV09IZ0CPDkGHKCIxZpGOcZvZNcDU8OFb7v78KVdqVuTuXcOvDfjoxHGtcpXAaqASuMfd/+8k7zkbmA2Qmpo6ae7cuRHHU1JSQkpys882AAAQHUlEQVRKSnSNaAFO1u6ScuftfZW8ubeCfUedNq1gSp9kzuuXzNBuSSQ14+U+9PtOLGp3ZKZPn77S3TMjKRtx4oiWmb0O9Knj0l3AX2omCjP7yN271fEe/dw918wGA4sILfG+vaG6MzMzfcWKFRHHmp2dTVZWVsTlW4ra7XZ33t95iCfe28OC9fspr6pmwoCuXJ/ZnyvG9SWlbcvYd1u/78SidkfGzCJOHA3tAFhMaPrtpy4RmqXbub7vdfcLT/K+B8wszd3zzCyN0JPodb1HbvjfHWaWDUwAGkwcEp0jZRU8uzKHJ97bw7b8Ejq3S+YLZw3ghin9yehT769YRBJUQ8uqd2qieucBXwXuCf/7Qu0CZtYNKHX342bWk9BtMu0B0ohyi6v56f+t47lVuZSWVzGuf1fuvXYsnx3bl/ZtNCtKROoW1L2He4C/m9mthLam/TyAmWUC33D324ARwENmVk1oBtc97r4xoHhbjOpq582tBfx56U6WbjtGm+QcZo3ry1fOOYOx6Z8aZhIR+ZRAEoe7FwIz6ji/Argt/PodYEyMQ2uxyiqqeGZlDo+8vZMdBUdJ7dyWa4a25l+uP58eKW2DDk9EmpGWMdop9froaDmPv7ubx9/dReHRcsamd+F3N4xn5pg03n5riZKGiERNiaOF2n+4jIff2sH/vr+H0vIqLsjozexpgzlrUHftnCcip0WJo4XJ+aiUP2Rv55kVOVS5M2tcX75x/pkM79NU8xxEJNEocbQQJxLG0yv2AnBdZn++Me1MPdktIo1OiaOZKyg+zgOLt/HEe7sBuH5yf27PGkLfru0DjkxEWioljmaquKyCOUt28OelOzleWc3nM9P51gVDlTBEpMkpcTQzlVXVPLl8L799bSuFR8v57Li+fP+iYQzq2THo0EQkQShxNCNvbi3gly9uYHvBUc4a1J1HLx+hh/ZEJOaUOJqBPYWl/NvLG3lt4wEG9ujAn76SyYUjemtarYgEQokjjpVXVvPgm9u5f/E2kpOMH146nFs/M4i2yVpHSkSCo8QRp5bvOsRPnlvHtvwSLh+bxs8uH0mfLu2CDktERIkj3pSWV3LPK5t5/N3d9Ovankdvmsz0jN5BhyUi8jEljjiyfNchfvD0GvYcKuXmqQO585LhdGijX5GIxBd9KsWBiqpqfvPqFuYs2UF6t/bM/drZnDW4R9BhiYjUSYkjYHsKS/nW3A9Ys7eIG6cM4KeXj6BjC9miVURaJn1CBWj+ujx+9MxaMPjDFycyc0xa0CGJiDQoKYhKzew6M9tgZtXhXf/qK3epmW0xs21m9uNYxtiUqqqde17ZzO1PrGJIagrzv32ekoaINBtB9TjWA1cDD9VXwMxaAQ8AFwE5wHIzm9fct48tKi3n23NXs2RrAV88awB3f3YUbZIDyd8iIqckqK1jNwENPfk8Bdjm7jvCZecCVwLNNnHsOniUmx59n31FZdxz9RhumDIg6JBERKJm7h5c5WbZwA/Ce43XvnYtcKm73xY+/jJwlrvfUc97zQZmA6Smpk6aO3duxHGUlJSQkpISfQOisKOoiv9ZVYY7fGdiO4Z2C/7p71i0Ox6p3YlF7Y7M9OnTV7p7vUMHNTVZj8PMXgf61HHpLnd/obHrc/c5wByAzMxMz8rKivh7s7OziaZ8tN7YdIB731hF704deOzmyQzuFR//ETd1u+OV2p1Y1O7G12SJw90vPM23yAX61zhOD59rVl5au4/vzF3NqL6d+fNXJ9OrU9ugQxIROS3xPCq7HBhqZoPMrA1wAzAv4Jii8uKaUNKYNKAb//u1s5U0RKRFCGo67lVmlgOcA7xsZgvD5/ua2XwAd68E7gAWApuAv7v7hiDiPRUvrM7lO3M/YNIZ3Xj05smk6KE+EWkhgppV9TzwfB3n9wEzaxzPB+bHMLRG8frGA3zvqdVMGdSdR26arPWmRKRFiedbVc3S6r1F3PHkKsb066KkISItkhJHI9p18Ci3Prac3p3a8WclDRFpoZQ4GslHR8u56dH3qXbnL7dMoWeKBsJFpGXSn8SNoLra+d7fV7OvqIwnZ5/NoJ4dgw5JRKTJqMfRCB5csp3sLQX87LMjmXRGt6DDERFpUkocp+n9nYf4r1e3csXYNL50ltaeEpGWT4njNBSWHOdbT65iQPcO/OfVYxpatFFEpEVQ4jgNv3hxIx8dreD+L0ygU7vWQYcjIhITShyn6M2tBby4Zh+3Tz+TUX27BB2OiEjMKHGcgrKKKn72f+sZ3LMj38w6M+hwRERiStNxT8HvF33InkOl/O/XzqJtcvD7aoiIxJJ6HFHaeqCYh97cwTUT0zn3zJ5BhyMiEnNKHFH61Sub6dg2mX+ZmRF0KCIigVDiiML63MO8sTmf2z4ziB5aUkREEpQSRxTuX7SNTu2S+erUgUGHIiISmKA2crrOzDaYWbWZ1bs5upntMrN1ZrbazFbEMsbatuwvZsGG/dx87kA665kNEUlgQc2qWg9cDTwUQdnp7n6wieNp0P2Lt9GxTStu+cygoEMREQlUUDsAbgKazRId2wtKeGntPr4+7Uy6dmgTdDgiIoGK9zEOB141s5VmNjuoIOa8uYO2yUncdp56GyIi5u5N88ZmrwN96rh0l7u/EC6TDfzA3escvzCzfu6ea2a9gdeAb7n7knrKzgZmA6Smpk6aO3duxLGWlJSQkpJS57XyKufbi0qZ3CeZW8e0rJlUJ2t3S6Z2Jxa1OzLTp09f6e71jjl/grsH9gVkA5kRlv0FoSTTYNlJkyZ5NBYvXlzvtZfX7vMzfvSSL/2wIKr3bA5O1u6WTO1OLGp3ZIAVHuFnd9zeqjKzjmbW6cRr4GJCg+ox9cLqXHp1asvZg3vEumoRkbgU1HTcq8wsBzgHeNnMFobP9zWz+eFiqcBSM1sDvA+87O4LYhnn4WMVLN5SwBVj02iV1DwG8kVEmlpQs6qeB56v4/w+YGb49Q5gXIxD+4SFG/ZTXlnNrHF9gwxDRCSuxO2tqngwb/U+BnTvwPj+XYMORUQkbihx1CO/uIx3th/kyvF9m83zJiIisaDEUY+X1+ZR7eg2lYhILUoc9Xh5bR4ZfToxNLVT0KGIiMQVJY46lFdWszbnMNOG9Qo6FBGRuKPEUYfN+49QXlXN2PQuQYciIhJ3lDjqsDbnMADj0jWbSkSkNiWOOqzNKaJbh9akd2sfdCgiInFHiaMOa3MOMza9q6bhiojUQYmjltLySrYeKGacxjdEROqkxFHLhn1HqHYYq/ENEZE6KXHUsmZvEQBj+6vHISJSFyWOWtbmHCatSzt6d2oXdCgiInFJiaOWtTlFen5DROQklDhqOFxawa7CUo1viIichBJHDWtzQ+MbevBPRKR+Qe0A+Gsz22xma83seTOr85PazC41sy1mts3MftzUcZ14YnyMblWJiNQrqB7Ha8Bodx8LbAV+UruAmbUCHgAuA0YCN5rZyKYMas3eIgb17EiX9q2bshoRkWYtkMTh7q+6e2X4cBmQXkexKcA2d9/h7uXAXODKpowr9MS4ehsiIidj7h5sAGYvAk+5+99qnb8WuNTdbwsffxk4y93vqOd9ZgOzAVJTUyfNnTs34hhKSkpo16Ejj20oZ2SPVpzbN5Ct2GOupKSElJSUoMOIObU7sajdkZk+ffpKd8+MpGyTfUKa2etAnzou3eXuL4TL3AVUAk+cbn3uPgeYA5CZmelZWVkRf292djZZWVlceMHpRtG8nGh3olG7E4va3fiaLHG4+4Unu25mNwFXADO87m5PLtC/xnF6+JyIiAQoqFlVlwI/BGa5e2k9xZYDQ81skJm1AW4A5sUqRhERqVtQs6ruBzoBr5nZajN7EMDM+prZfIDw4PkdwEJgE/B3d98QULwiIhIWyCiwuw+p5/w+YGaN4/nA/FjFJSIiDdOT4yIiEhUlDhERiYoSh4iIREWJQ0REohL4k+NNwcwKgN1RfEtP4GAThRPP1O7EonYnlmjbfYa794qkYItMHNEysxWRPmrfkqjdiUXtTixN2W7dqhIRkagocYiISFSUOELmBB1AQNTuxKJ2J5Yma7fGOEREJCrqcYiISFSUOEREJCoJkzjM7FIz22Jm28zsx3Vcb2tmT4Wvv2dmA2MfZeOLoN3fN7ONZrbWzN4wszOCiLMpNNT2GuWuMTM3sxYxZTOSdpvZ58O/9w1m9r+xjrEpRPDf+gAzW2xmH4T/e59Z1/s0J2b2iJnlm9n6eq6bmd0X/pmsNbOJjVKxu7f4L6AVsB0YDLQB1gAja5W5HXgw/PoGQtvZBh57DNo9HegQfv3NltDuSNseLtcJWAIsAzKDjjtGv/OhwAdAt/Bx76DjjlG75wDfDL8eCewKOu5GaPc0YCKwvp7rM4FXAAPOBt5rjHoTpccxBdjm7jvcvRyYC1xZq8yVwF/Cr58BZpiZxTDGptBgu919sf9jM61lhHZabAki+Z0D/BvwK6AslsE1oUja/TXgAXf/CMDd82McY1OIpN0OdA6/7gLsi2F8TcLdlwCHTlLkSuBxD1kGdDWztNOtN1ESRz9gb43jnPC5Ost4aBOpw0CPmETXdCJpd023EvrrpCVosO3hbnt/d385loE1sUh+58OAYWb2tpktC+/I2dxF0u5fAF8ysxxC+/x8KzahBSraz4CIBLKRk8QfM/sSkAmcH3QssWBmScB/AzcFHEoQkgndrsoi1MNcYmZj3L0o0Kia3o3AY+7+X2Z2DvBXMxvt7tVBB9bcJEqPIxfoX+M4PXyuzjJmlkyoK1sYk+iaTiTtxswuBO4itAf88RjF1tQaansnYDSQbWa7CN3/ndcCBsgj+Z3nAPPcvcLddwJbCSWS5iySdt8K/B3A3d8F2hFaCLAli+gzIFqJkjiWA0PNbJCZtSE0+D2vVpl5wFfDr68FFnl4dKkZa7DdZjYBeIhQ0mgJ97pPOGnb3f2wu/d094HuPpDQ+M4sd18RTLiNJpL/1v+PUG8DM+tJ6NbVjlgG2QQiafceYAaAmY0glDgKYhpl7M0DvhKeXXU2cNjd8073TRPiVpW7V5rZHcBCQrMvHnH3DWb2S2CFu88D/kyo67qN0GDTDcFF3DgibPevgRTg6fBcgD3uPiuwoBtJhG1vcSJs90LgYjPbCFQBd7p7s+5dR9jufwb+ZGbfIzRQflNz/+PQzJ4k9EdAz/DYzd1AawB3f5DQWM5MYBtQCtzcKPU285+biIjEWKLcqhIRkUaixCEiIlFR4hARkagocYiISFSUOEREJCpKHCIiEhUlDhERiYoSh0gMmNlAM9tsZk+Y2SYze8bMOgQdl8ipUOIQiZ3hwB/cfQRwhNAeMCLNjhKHSOzsdfe3w6//BnwmyGBETpUSh0js1F7fR+v9SLOkxCESOwPC+0AAfAFYGmQwIqdKiUMkdrYA/2Rmm4BuwB8DjkfklCTEsuoicaLS3b8UdBAip0s9DhERiYr24xARkaioxyEiIlFR4hARkagocYiISFSUOEREJCpKHCIiEpX/Dy/HENlO0aXLAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(p, np.log10(odds),label=\"log10odds\")\n",
    "plt.xlabel(\"p\")\n",
    "plt.ylabel(\"log10odds\")\n",
    "plt.legend()\n",
    "plt.title(\"log10odds\")\n",
    "plt.grid()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "底数不是特别重要，但是合适的底方便对比。"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2019-02-25T01:47:42.673659Z",
     "start_time": "2019-02-25T01:47:42.415586Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEWCAYAAABxMXBSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xl8VNX5+PHPkxCWhJ2EEEjYJGyyBAgIIhhALCJiUdyXUq30q/JtqdrWttra1v5qtVZrXSr161oVcKugCLINiCCbArIk7JAFSAIESCAJSZ7fHzPQGBIyIZm5k5nn/Xrlxdx7T+55zkzIk3PPveeIqmKMMcZ4K8zpAIwxxtQvljiMMcbUiCUOY4wxNWKJwxhjTI1Y4jDGGFMjljiMMcbUiCUOE5JEZK+IXOF0HN4SERWRblUcmyIiK/wdkwldljiMuUAiMkpElorIMRHZW8nxzp7jJ0UktT4lKmPOxxKHMReuAHgV+HkVx98FvgHaAL8B3heRGD/FZozPWOIwIU1EGonIsyKS5fl6VkQalTv+CxE54Dn2o/KXjFR1jaq+Beyu5LzdgYHA71T1lKp+AHwLXO9lvT8vV+9dFc7dRkTmiMhxEVkDXFTumIjIMyKS7Tn+rYj0qdt3zYQ6Sxwm1P0GGAokAf2BIcAjACIyDngAuALoBqTU4LwXA7tV9US5fRs9+72p9yFgLJDoqb+8F4BCIA64y/N1xpXASKA70AK4EThcg7iNqZYlDhPqbgP+oKrZqpoD/B64w3PsRuA1Vd2iqieBx2pw3qbAsQr7jgHNalDvZlUtKF+viITj7rX8VlULVHUz8Ea5Ok576ugJiKpuU9UDNYjbmGpZ4jChrj2wr9z2Ps++M8fSyx0r/7o6+UDzCvuaA2d6IDWpt3y5GKBBVcdVdQnwPO5eSbaIzBCRinEYUyuWOEyoywI6ldvu6NkHcACIL3csoQbn3QJ0FZFm5fb19+z3pt6ECsfOyAFKznMcVX1OVQcBvXFfsqpq8N6YC2KJw4S6d4FHRCRGRKKB3wL/9hybDfxQRHqJSCTwaPlvFJEwEWkMRLg3pbGINARQ1e3ABuB3nv2TgH7AB17WO0VEenvq/d2ZOlW1FPgQeExEIkWkN/CDcjENFpFLRCQC911fhUBZnbxTxnhY4jCh7nFgHbAJ911PX3v2oaqfAc8BS4GdwFee7yny/DsSOAXMw/1X/yng83LnvhlIBo4CTwCTPeMZ3tT7LLDEU++SCjFPwz2GchB4HXit3LHmwL88de7DPTD+VA3eD2OqJbaQkzHeEZFewGagkaqWOB2PMU6xHocx5yEikzzPXLQC/gLMtaRhQp0lDmPO78dANrALKAXudTYcY5xnl6qMMcbUiPU4jDHG1EgDpwPwhejoaO3cubPX5QsKCoiKivJdQAHK2h1arN2hpabtXr9+fa6qejUJZ1Amjs6dO7Nu3Tqvy7tcLlJSUnwXUICydocWa3doqWm7RWRf9aXc7FKVMcaYGrHEYYwxpkYscRhjjKmRoBzjqMzp06fJyMigsLDwnGMtWrRg27ZtDkTlrLpod+PGjYmPjyciIqKOojLGBLqQSRwZGRk0a9aMzp07IyLfOXbixAmaNWtWxXcGr9q2W1U5fPgwGRkZdOnSpQ4jM8YEMscuVYlIgogsFZGtIrJFRH5aSRkRkedEZKeIbBKRgRdaX2FhIW3atDknaZgLJyK0adOm0l6cMSZ4OdnjKAEeVNWvPWsWrBeRhaq6tVyZq3AvnZkIXAK85Pn3gljSqHv2nhoTehzrcajqAVX92vP6BLAN6FCh2LXAm+r2FdBSROL8HKoxxgS8IwXFnCou9UtdATHGISKdgQHA6gqHOvDdJTIzPPvOWUNZRKYCUwFiY2NxuVzfOd6iRQtOnDhR8dsAKC0trfJYXYqLi+PAgQtb/vnll1/mxRdfZM+ePezZs4c2bdoA7nGGX/ziF3z++edERkby0ksvkZSU5NU5S0tLue222xg3bhzf//73v3Psiy++4LnnnuO9996r9jyFhYXnvN+BLD8/v17FW1es3cHtra1FrMoq4fkxkYSJ+LTdjicOEWmKe1W06ap6/ELPo6ozgBkAycnJWvGJyW3btlU5EOzPwfELrWfMmDFMnjyZlJQUmjZtevY88+bNY9++fezatYvVq1fz05/+lNWrK+bfyp04cYKIiAiaNGlyTlyRkZE0aNDAq3gbN27MgAEDat4ohzj1JPGi3YsAuKLrFX6vG+wJ6mD3zJYv6ZsQxuhRwwDfttvR5zg8y1t+ALytqh9WUiST766tHO/ZV6+pKj//+c/p06cPffv2ZdasWQCUlZVx33330bNnT8aOHcv48eN5//33ARgwYACVzb/18ccfc+eddyIiDB06lLy8PA4cOFBlHarKtGnT6NGjBxMnTiQ7O/vsuebPn0/Pnj0ZOHAgH374349j2bJlJCUlkZSUxIABA/zSOwtGjy9/nMeXP+50GCYIFZeUsS3rOP0TWvqlPsd6HOIeVf0/YJuq/q2KYnOAaSIyE/eg+DFVvbBrPeX8fu4Wtmb9t3NTWlpKeHh4rc7Zu31zfnfNxV6V/fDDD9mwYQMbN24kNzeXwYMHM3LkSL788kv27t3L1q1byc7OplevXtx1113nPVdmZiYJCf/NrfHx8WRmZrJy5cpK61i1ahVpaWls3bqVXbt2MWTIEO666y4KCwu55557WLJkCd26deOmm246e86//vWvvPDCCwwfPpz8/HwaN258YW+SMcYn0g6eoLi0jH7xLfxSn5M9juHAHcBoEdng+RovIv8jIv/jKTMP2I173eV/Afc5FGudWrFiBbfccgvh4eHExsZy+eWXs3btWlasWMENN9xAWFgY7dq1Y9SoUXVex/Lly8/uj4uLY/To0QCkpqbSpUsXEhMTERFuv/32s+caPnw4DzzwAM899xx5eXk0aOD4FU5jTDkbMvIA6B8f5D0OVV0BnPdeTnWvMnV/XdddsWdQnx8A7NChA+np/71/ICMjgw4dKt6cVjsPP/wwV199NfPmzWP48OEsWLCAnj171mkdxpgLtyk9j9ZRDYlv1cQv9dlcVQ4YMWIEs2bNorS0lJycHJYvX86QIUMYPnw4H3zwAWVlZRw6dMirOyImTpzIm2++iary1Vdf0aJFC+Li4qqsY+TIkWf3Hzx4kKVLlwLQs2dP9u7dy65duwB49913z9axa9cu+vbtyy9/+UsGDx5MamqqT94XY8yF2ZRxjH7xLfz2XJVdc3DApEmTWLVqFf3790dEePLJJ2nXrh3XX389ixcvpnfv3iQkJDBw4EBatHBfs3zuued48sknOXjwIP369WP8+PG88sorjB8/nnnz5tGtWzciIyN57bXXzlvHpEmTWLJkCb1796ZDhw4MG+a+A6Nx48bMmDGDq6++msjISEaMGHF2EPzZZ59l6dKlhIWFcfHFF3PVVVc588bVcy9PeNnpEEwQKigqYUf2Cb7Xp53/KlXVoPsaNGiQVrR169Zz9p1x/PjxKo/524kTJ1RVNTc3V7t27aoHDhzwWV111e7zvbeBaOnSpU6H4Ahrd3D6aleudvrlJ7po68Hv7K9pu4F16uXvWOtxBJgJEyaQl5dHcXExjz76KO3a+fGvCONTc9PmAnBNj2scjsQEk00ZxwDo56eBcbBLVQEnFJ5wDVVPr3oasMRh6tbGjDw6tGxCTLNGfqszpAbH3b0xU5fsPTXGWWcGxv0pZBJH48aNOXz4sP2iq0PqWY/DHgg0xhlHC4rZf+SkXy9TQQhdqoqPjycjI4OcnJxzjhUWFobkL7+6aPeZFQCNMf638eyDf/7tcYRM4oiIiKhylTqXy1WvJumrK6HabmOCxZmB8T6WOIwJTm9NesvpEEyQ+Wb/UbrGRNG8cYRf67XEYYyfJLRIqL6QMV46dvI0X+48zG1DO/q97pAZHDfGabM2z2LW5llOh2GCxPwtByguLeP7SXU7N503rMdhjJ+8tO4lAG7qc1M1JY2p3scbsujcJtLvt+KC9TiMMabeOXS8kFW7DzMxqYPfJjYszxKHMcbUM3M3ZqEKE/u3d6R+SxzGGFPPzNmYRZ8OzenWtqkj9Tu95virIpItIpurOJ4iIsfKrRD4W3/HaIwxgWR3Tj6bMo5xbX//D4qf4fTg+OvA88Cb5ynzhapO8E84xvjO+ze+73QIJgjM2ZiFCEzoH+dYDI4mDlVdLiKdnYzBGH+Jjox2OgRTz5WWKR9+ncmQzq2Ja+GfZWIrI05P+udJHJ+oap9KjqUAHwAZQBbwkKpuqeI8U4GpALGxsYNmzpzpdQz5+fk0berMtUInWbv9a/7B+QCMazfO73WDfd7BYFVWCS9vKuL+pEYMbnf+v/tr2u5Ro0atV9Vkrwp7u+KTr76AzsDmKo41B5p6Xo8HdnhzzspWADyfYF8hrCrWbv+6/LXL9fLXLnekblX7vOu70tIyHfO0S6/82zItLS2rtrwvVwAM6LuqVPW4quZ7Xs8DIkTE+vvGmJDz2eaD7MzOZ9roboSF+f/ZjfICOnGISDvxPN0iIkNwx3vY2aiMMca/ysqUfyzZQdeYKMb3dW5Q/AxHB8dF5F0gBYgWkQzgd0AEgKr+E5gM3CsiJcAp4GZPl8oYY0LGwm2HSD14gmdu6k+4w70NcP6uqluqOf487tt1jTEmJJWVKc8t3kGnNpFc08+ZJ8Urcvo5DmNCxrzb5jkdgqmHZq5NZ0vWcf52Y38ahAfG6IIlDmP8JDIi0ukQTD2TfaKQJz7bxiVdWjNpgHNPilcUGOnLmBDw4toXeXHti06HYeqRxz/ZRuHpMv40qa8js+BWxRKHMX4ye8tsZm+Z7XQYpp5Ytj2HORuzuDflIscmM6yKJQ5jjAkwJ4tLePQ/m+kaHcW9KRc5Hc45bIzDGGMCiKryyEebST96knfvGUrjiHCnQzqH9TiMMSaAzF6XzoffZDJ9THeGdm3jdDiVssRhjDEBYmvWcX778RZGJEYzbXQ3p8Opkl2qMsZPXFNcTodgAtjxwtPc/87XtIyM4JmbkgLiCfGqWOIwxhiHFZeUce+/15N+5CTv3DOU6KaNnA7pvCxxGOMnf135VwAeuvQhhyMxgURV+eUHm/hy52GevqE/Q7q0djqkatkYhzF+8sn2T/hk+ydOh2ECzJML0vjom0weurI71w+Kdzocr1jiMMYYh7zyxW5ecu3i1ks6cv+owB0Mr8guVRljjANeXraLP3+Wyvi+7fjDxIsDakqR6ljiMMYYP3vRtZMn56cxoV8cz96UFDCz3nrLEocxftIkoonTIRiHqSrPLNrBc4t3cG1Se56+IXCmSq8Jp1cAfBWYAGSrap9Kjgvwd2A8cBKYoqpf+zdKY+rGZ7d95nQIxkGnS8v49Yff8t76DG4YFM8T1/cL6Gc1zsfpVPc6MO48x68CEj1fU4GX/BCTMcbUqfyiEu5+Yx3vrc/gp2MSeXJy/U0a4PzSsctFpPN5ilwLvOlZZ/wrEWkpInGqesAvARpTh/647I8APHr5ow5HYvxp/+GTTH1rHTuy8/nL9X25aXBHp0OqtUAf4+gApJfbzvDsOydxiMhU3L0SYmNjcblcXleSn59fo/LBwtrtXx9s+ACAETrC73WDfd5O+DanhH9uKgJg+oBGxBbsxuXa7Ze6fdnuQE8cXlPVGcAMgOTkZE1JSfH6e10uFzUpHyys3f7Vcm9LAMfec/u8/aesTHlp2S7+9nUaPWKbMeOOZDq28e/Swb5sd6Anjkwgodx2vGefMcYEpOzjhTwweyMrduZyTf/2/OX6vkQ2DPRftTUT6K2ZA0wTkZnAJcAxG98wxgSqpWnZPDR7IwXFJfz5ur7cPDihXj3Y5y2nb8d9F0gBokUkA/gdEAGgqv8E5uG+FXcn7ttxf+hMpMbUXpvIwFyUx9ReflEJf/p0G++u2U+P2GbMvHUoibHNnA7LZ5y+q+qWao4rcL+fwjHGpz648QOnQzA+sHJXLj9/bxNZx04xdWRXHhjbPSCXe61LgX6pyhhjAlLeyWL+PC+VWevS6RIdxfv/M4xBnQJ/SvS6YInDGD/51aJfAfDnK/7scCSmNlSV/2zI5PFPtpF36jQ/vrwr08d0p0nD4O5llGeJwxg/WZWxyukQTC1tzTrOY3O3sGbPEfontOStSX3p3b6502H5nSUOY4ypxtGCYv62cDtvr95HiyYRPP79PtwypGO9njakNixxGGNMFQpPl/LGyr08v3QnBUUl3DG0Ez8b252WkQ2dDs1RljiMMaaC0jLl4w2ZPP35djLzTjGqRwy/Gt+L7kF8i21NWOIwxk/im9eP9aRDmaqyYMshnv48jR3Z+fTp0JynJvfj0m7RTocWUCxxGOMn/77u306HYKqgqizels3fF+/g28xjdI2J4oVbB3JVn3aEheg4xvlY4jDGhKyyMmXRtkM8t2QHmzOP07F1JE9O7sd1AzrUy5X5/MUShzF+Mn3+dACeHfesw5GYktIyPtl0gBddO9l+KJ9ObSJ5anI/vj+gAxGWMKplicMYP9lwcIPTIYS8gqIS3luXzisr9pBx9BQ9Ypvx95uTuLpvnPUwasAShzEm6B06Xsibq/by76/2c+zUaZI7teK3E3pzRa9YG8O4AJY4jDFBa0N6Hv/cWMi6z5dQqsq4i9vxoxFdGdSpldOh1WuWOIwxQaWopJR53x7gjZX72JCeR+Nw+MGlXfjBsM5+X4UvWFniMMZPurfp7nQIQS39yEneWbOf2WvTOVxQTNeYKB67pjexp/Zy1RW9nQ4vqDi9kNM44O9AOPCKqj5R4fgU4Cn+u1zs86r6il+DNKaOzLhmhtMhBJ3TpWUs3pbNu2v2s3xHDgKM6RXLD4Z1Zni3NogILtc+p8MMOo4lDhEJB14AxgIZwFoRmaOqWysUnaWq0/weoDEmYO3JLWDW2nTeX59Bbn4Rsc0b8ZPRidw8JIG4Fk2cDi/oOdnjGALsVNXdAJ51xa8FKiYOY4LC1LlTAet5XKiCohI+/fYA76/LYM3eI4QJjO7ZlpsHdySlR4zdTutH4l6d1YGKRSYD41T1R57tO4BLyvcuPJeq/gzkANuBn6lqehXnmwpMBYiNjR00c+ZMr2PJz8+nadOmF9iS+sva7V/TN3geAExy5gHA+vh5l6mSdqSMFZklrDtUQlEpxEYKI+IbMLx9A1o1rj5Z1Md214WatnvUqFHrVTXZm7KBPjg+F3hXVYtE5MfAG8Doygqq6gxgBkBycrKmpKR4XYnL5aIm5YOFtdu/Wu5tCeDYe16fPu+d2Sf48OtMPt6QRWZeIc0aNWDSwASuHxRPcqdWiHj/7EV9andd8mW7vUocIhIFnFLVMhHpDvQEPlPV07WoOxNIKLcdz38HwQFQ1cPlNl8BnqxFfcaYAHbwWCGfbMriPxsy2Zx5nPAwYURiNL8Y14PvXdyOxhGhszRroPO2x7EcGCEirYDPgbXATcBttah7LZAoIl1wJ4ybgVvLFxCROFU94NmcCGyrRX3GmABztKCYzzYfZO7GLL7acxhV6Bffgt9O6M01/dsT06yR0yGaSnibOERVT4rI3cCLqvqkiNRq4h1VLRGRacAC3LfjvqqqW0TkD8A6VZ0D/EREJgIlwBFgSm3qNMZJSe2SnA4hIBwvPM3nWw7xyaYsVuzIpaRM6RodxU/HJDKxf3u6xoTeeER943XiEJFhuHsYd3v21brfqKrzgHkV9v223OtfAb+qbT3GBIJQnhX3eOFpFm87xKebDrB8ey7FpWV0aNmEu0d04Zp+7bm4ffMajVsYZ3mbOKbj/gX+kadX0BVY6ruwjDH1Xd7JYhZty+azbw/wxQ53sohr0Zg7h3Xi6n5xJCW0tGRRT3mVOFR1GbCs3PZu4Ce+CsqYYHT7h7cDwb0SYPaJQhZuPcT8zQdZteswJWVKh5ZNuHNYJ8b3iyMpvqXNRhsEzps4RGQuUOWDHqo6sc4jMiZIZRzPcDoEn9h3uIDPtxxiwZaDrN9/FFXo3CaSH43oylV92tEvvoX1LIJMdT2Ov3r+vQ5oB5z5U+kW4JCvgjLGBC5VZXPmcRZuPcjnWw+RevAEAD3bNWP6mO6M69OO7rFNLVkEsfMmDs8lKkTk6QpPFM4VkXU+jcwYEzCKSkr5avcRFm09xKJthzhwrJAwgcGdW/PohN5c2TuWhNY2ZXmo8HZwPEpEupabV6oLEOW7sIwxTjtSUMzS1GwWpx5iWVoOBcWlNI4IY2RiDA+M7c6YXrG0jmrodJjGAd4mjp8BLhHZDQjQCfixz6IyJggNix/mdAjnpapsP5TP4tRDLNmWzdf7j1Km0LZZIyYmdWBs77ZcelG0PcFtvL6rar6IJOKeagQgVVWLfBeWMcHnz1f82ekQzlF4upRVuw+zZFs2S1Kzycw7BcDF7ZszbXQiV/RqS5/2LexOKPMd1d1VdV0Vhy4SEVT1Qx/EZIzxoYyjJ1maloMrNZsvd+VSeLqMJhHhDO8WzbTR3RjVoy3tWjR2OkwTwKrrcVzj+bctcCmwGPelqlHASsAShzFeun729QB8cOMHfq33dGkZ6/YeZVZaMf/vm2VsP5QPQELrJtyUnMDoXrFc0qW1XYIyXqvurqofAojI50DvMxMOikgc8LrPozMmiBw+ebj6QnXkwLFTLEvLwZWWw5c7czlRVEK4wNCLmnNjcgKjerala3SU3TJrLoi3g+MJ5WapBfczHB19EI8x5gIUl5Sxbt+Rs8ki7ZD72Yq4Fo2Z0D+OlB5tKTuwjauuGOpwpCYYeLvW4mIRWSAiUzyr8n0KLPJdWMYEr1OnTnH55ZdTWloKwBtvvEFiYiKJiYm88cYblX7PY489RocOHUhKSiIpKYl58+aRfuQk//5qHz96Yx0D/vA5t/5rNa9+uYc2TRvyq6t68ki/Io68fh/v/eI6vpn7Ok0auHsXN998Mzt27PBbe03w8fauqmkiMgkY6dk1Q1U/8l1YxgSvV199leuuu47w8HCOHDnC73//e9atW4eIMGjQICZOnEirVq3O+b5p//tThk6awvLtOfx1aw67l7vnGY1v1YRJAzswMjGGS7tF07RRA0pLS+ne/UoWLlxIfHw8gwcPpl27dgDce++9PPnkk/zrX//ya7tN8KjJ0rErca+LocAa34RjTPAa02UMAG//6m3eeecdABYsWMDYsWNp3bo1AGPHjmX+/PnccsstZ5+rWL49hw+/zuDASYjKW0vjiDCGdm3D7UM7MbJ7DBfFnDtWsWbNGrp160bXrl0Bdy/jyy+/ZMqUKYwYMYIpU6ZQUlJCgwaBvnq0CUTeLh17I/AU4MJ9V9U/ROTnqvq+D2MzJqg8evmjFBcX88LuF+jcuTMAmZmZJCT8dwXlNm3jWLh2G1832MgXO3I5eLwQgLCiEkq+XUh4+kouGXYJzzzwt0p7JWdUPG98fDxr1rj/3gsLC6Nbt25s3LiRQYMG+aClJth5++fGb4DBqpoNICIxuMc4apU4RGQc8Hfci0K9oqpPVDjeCHgTGAQcBm5S1b21qdMYJ+Xm5tKyZcuz22VlSubRk/zt8zSW78hl2bKdEN6IhOYHuSwxmpGJMYzsHkN4UTLR0f+HiPDoo4/y4IMP8uqrr15wHG3btiUrK8sSh7kg3iaOsDNJw+Mw3g+sV0pEwoEXgLFABrBWROao6tZyxe4GjqpqNxG5GfgL7rXOjal3rnr7Kk4XnKbg5CneXr2PL7bn8unqw+Tt2khMw+EkJbSkZ7MSJnwvhV/eP5YG4eX/izU5++qee+5hwoQJ562rQ4cOpKenn93OyMggOjr67HZhYSFNmjSp7FuNqZa3iWO+iCwA3vVs30SFJV8vwBBgZ7mJE2cC1wLlE8e1wGOe1+8Dz4uIqGqVa4QYU53pG6bTcm/L7+y78eIbuW/wfZw8fZLxb48/53umJE1hStIUck/mMnn25HOO35t8Lzf1uYn0Y+nc8dEd3zlWVhbBVR1/Qkb6pZws6MiBvIf49Xtf06BxIVG9G5DtWsEfr72H8T160+eJL5h3VToL3/ruYk8P9HmAickTWZm+kjt/dSfHo46T8noKRUeLSP1XKkuXLCWpXRKLdi/i8eWPo6XKmm/WcMlTl9CoVSOy387mlw/+8uz5tm/fTp8+fWrzNpoQ5u1dVT/3TD9ymWdXXdxV1QFIL7edAVxSVRlVLRGRY0AbILfiyURkKjAVIDY2FpfL5XUg+fn5NSofLEK13aWlpeTl5X1n3/bt23EVuCgsLTznGEBqaiquPBfHTh+r9PiWrVtw5brILszm6NE8SkviKCnqRklxIiXFnXh5exjQh7CIXTRNjCbs+B9p2qEMEWh7ZUt+OmkqvwhrxJU3XMk3Jd9AHqTPTKfNpW2I7BjJ4w89zvT90ykqKyK/aT7xN8STl5fHyfSTlJSVsG7dOvKa5rHx6Maz8cV9P44NT22AMrhm/DXExMTgcrk4cuQIJSUlpKamkpqa6pP3OJCE6s+5L9st3vzxLiKxuH+JA2Sqaq0XcRKRycA4Vf2RZ/sO4BJVnVauzGZPmQzP9i5PmXMSR3nJycm6bp33y4W4XC5SUlJq3oh6ztpdd3Lzi1ixI5fl23NYviOX3Hz3HKA92zVjZPcYRibG8OsvbiIsrJS/9fsbzzzzDG+99Vat633++efp2LEjEydWvxjnmXY/88wzNG/enLvvvrvW9dcH9nPuHRFZX2HdpSpVN8lhEvBPoAWQ6dkdLyJ5wH2q+rXXUZ0rE0gotx1fro6KZTJEpIEnDv/N22BMFYpLyvh6/1FPoshhc+ZxAFpHNeSybtGM7B7DiMRoYpv/d7LAsC/dD/wNHDiQUaNGUVpaSnh47eaHmjZtWvWFKmjZsiV33HFH9QWNqUJ1l6peB36sqqvL7xSRocBrQP9a1L0WSPQsCpUJ3AzcWqHMHOAHwCpgMrDExjeMU/YfPsmyHTksS8th1a5cCopLaRAmDOzYioeu7M7I7jHnnYJ8Qvf/Dmjfdddd/gr7HD/84Q8dq9sEh+oSR1TFpAGgql+JSK1WAPSMWUwDFuC+HfdVVd0iIn8A1qnqHOD/gLdEZCdwBHdyMcYvThaXsHr3EZZtz2HZ9hz25BYA7ie1vz8BHa+AAAAW7UlEQVSgAyO7x3DpRW1o1jjCq/M9dOlDvgzXGL+pLnF8JiKf4n6W4sxAdgJwJzC/tpWr6jwq3J2lqr8t97oQuKG29RjjDVVlZ3b+2USxes8RikvKaBwRxrCubfjBMPeT2l1sVlkT4qqbVv0nInIV7ttizw6OAy94fukbU6/lF5Xw5c5cXGk5LN+ec3YFvMS2TbljaCdSesQwuHPdrFWR8noKAK4prlqfyxgnVXs7rqp+Bnzmh1iM8TlVJeNEGS8v24UrLYe1e49QUqY0bdSA4d3acP+oblzeI4YOLe3hOGOqcsEznInIDFWdWpfBGOMLBUUlrNx1mKVp2bhSs8k6Vgik0rNdM+4e0YWU7m0Z1KkVDRvUajIEY0JGdbfjtq7qEHDu47XGBIi9uQUsTctmSWo2q3cfobi0jKiG4VyWGM334suYOvEy4lpYr8KYC1FdjyMH2Ic7UZyhnu22vgrKmJo6XVrG2r1HWLItmyVp2ezOcd8B1TUmijuHdWJ0z7Ykd25NwwZhuFwuSxrG1EJ1iWM3MEZV91c8ICLplZQ3xm+OFhTj2p7Nom3ZLE/L4URRCQ3Dw7ika2vuHNqJ0T1j6dgm0ukwz7rx4hudDsGYOlFd4ngWaAWckziAJ+s+HGPOb1dOPou2HmLxtmzW7TtCmUJMs0aM7xvH6F5tuaxbNFGNAnNxovsG3+d0CMbUiepux33hPMf+UffhGPNdpWXK1/uPsnDrIRZtPcRuz0N4veOaM21UN8b0iqVvh6qf1g4kJ0+fBCAyInB6QcZcCG9XALyukt3HgG8rrNNhTK0Vni5lxY5cPt96kMXbsjlcUExEuDC0axt+OLwzo3vF1svbZc9M127PcZj6zts+/d3AMGCpZzsFWA90EZE/qGrtp/k0Ie144WmWpmYzf/NBlm3P4WRxKc0aNSClZ1uu7B3L5T1iaO7l1B7GGN/yNnE0AHqdmU7dM836m7jXz1gOWOIwNXY4v4iFWw8xf8tBvtyZy+lSJaZZIyYN6MD3Lm7H0K5t7NkKYwKQt4kjocIaHNmefUdE5LQP4jJBKvtEIQs2H2TetwdZvecwZQodW0fyw+Fd+N7F7RiQ0LJejFcYE8q8TRwuEfkEeM+zPdmzLwo4dzk0Y8rJOVHE/M0H+GTTAdbsPYIqXBQTxf2jujGuTzt6xzW3SQONqUe8TRz3A+WXjn0D+MCzNsYoXwRm6rejBcV8tvkgczdmne1ZdGvblJ+MTuTqfnEktm0acsliStIUp0Mwpk54u+a4isgKoBj3k+NrbEElU1F+UQkLtx5kzoYsvtiRS0mZ0iU6immjunF1v/b0aNfM6RAdZYnDBAtvb8e9EXgKcOGebuQfIvJzVX3fh7GZeuB0aRnLt+fwnw1ZLNx6kMLTZXRo2YS7R3Thmn7tubi9XYY6I/dkLgDRkdEOR2JM7Xh7qeo3wOAzz2yISAywCLigxOGZPHEW0BnYC9yoqkcrKVcKfOvZ3K+qEy+kPlO3VJWNGcf46OsM5m46wJGCYlpFRnDDoASuTWrPwI6tbIC7EpNnTwbsOQ5T/3mbOMIqPOh3GKjNfZIPA4tV9QkRediz/ctKyp1S1aRa1GPq0MFjhXz4TQYfrM9gV04BDRuEMbZ3LNd5llGNCLdbZ40JBd4mjvkisgB417N9ExWWfK2ha3E/RAjugXYXlScO47CiklIWbc1m9rp0vtiRQ5lCcqdW/Pm6rlzdL84eyjMmBIm3Y9wicj0w3LP5hap+dMGViuSpakvPawGOntmuUK4E2ACUAE+o6n/Oc86pwFSA2NjYQTNnzvQ6nvz8fJo2bVqzRgSB87U7M7+MZemnWZlVQv5paN1YGN6+AZd1aEBsVP3uWTj1eU/fMB2AZ5Oe9XvdYD/noaam7R41atR6VU32pqzX04iq6gfAB96WF5FFQLtKDv2mwnlVRKrKXp1UNVNEugJLRORbVd1VRXwzgBkAycnJmpKS4m2ouFwualI+WFRsd+HpUuZ9e4C3V+9n/b6jRIQLV/aO48bBCVzWLZrwIBm3cOrzbrnX/beRUz9r9nMeWnzZ7upWADyB+/bbcw7h/p3fvKrvVdUrznPeQyISp6oHRCQO95PolZ0j0/PvbhFxAQOAShOHuXB7cwv491f7eP/rDPJOnqZLdBS/Ht+T6wfG06ZpI6fDCxr3Jt/rdAjG1InqplX31Y33c4AfAE94/v24YgERaQWcVNUiEYnGfZnM1gCpI2VlysacEl5/bQ2utBwahAnfu7gdt13SkWEXtbFbaH3gpj43OR2CMXXCqRVvngBmi8jduJemvRFARJKB/1HVHwG9gJdFpAz3HVxPqOpWh+INGgVFJby/PoPXvtzD3sNFxDQ7zvQrErl1SEfaNm/sdHhBLf2Ye9HMhBYJDkdiTO04kjhU9TAwppL964AfeV6vBPr6ObSgdeDYKV7/ci/vrNnPicISBnRsybj4Uh64YbTNQOsnd3x0B2DPcZj6LzDX2DR1JvXgcWYs282cjVkocFWfdtx9WRcGdGyFy+WypGGMqTFLHEFq/b6jvLh0J4tTs4lsGM6dwzrzw+GdSWhty5YaY2rHEkeQWbXrMM8t3sGq3YdpFRnBA2O7c+ewTrSMbOh0aMaYIGGJI0h8tfswzyzczuo9R4hp1ohHru7FrZd0JLKhfcTGmLplv1XquY3pefz18zS+2JFL22aN+N01vbllSEcaR4Q7HZqp4MFhDzodgjF1whJHPbU7J58n56cxf8tBWkc15JGre3H70E6WMALYNT2ucToEY+qEJY565nB+EX9fvIN3Vu+nUYMwfnZFd+4e0YWmjeyjDHRpuWkA9Iju4XAkxtSO/bapJ4pLynhj5V7+vngHp06XcsuQBH46pjsxzWxKkPrix5/8GLDnOEz9Z4mjHli+PYfH5m5hd04BKT1ieOTq3nRrG3qzfRpjAoMljgB26Hghf5i7lU+/PUDnNpG8OiWZ0T1jnQ7LGBPiLHEEoLIy5e3V+3hyfhpFpWU8OLY7Uy/vSqMGNvBtjHGeJY4Asze3gF+8v4k1e48wvFsbHv9+X7pERzkdljHGnGWJI0CUlSlvfbWPJz5LpUG48NTkfkweFG/TmweRR0Y+4nQIxtQJSxwBIPtEIQ/O3sgXO3IZ2T2Gv1zfl7gWTZwOy9SxK7pWubaZMfWKJQ6HLduew4OzN5BfVMKfJvXh1iEdrZcRpDYc3ABAUrskhyMxpnYcmVNbRG4QkS0iUuZZvKmqcuNEJE1EdorIw/6M0ddKSst44rNUfvDqGtpENWLutMu47ZJOljSC2PT505k+f7rTYRhTa071ODYD1wEvV1VARMKBF4CxQAawVkTmBMMqgEcLivnfd79hxc5cbr2kI7+d0NumCjHG1BtOrQC4Dajur+shwE5V3e0pOxO4FqjXiWNL1jF+/NZ6so8X8eT1/bhxsC0jaoypXwJ5jKMDkF5uOwO4pKrCIjIVmAoQGxuLy+XyuqL8/Pwalb9QG7JLeHFjEVENhIcHN6JtwS5crl0+r7cq/mp3oHGq3Xl5eQCOvef2eYcWX7bbZ4lDRBYB7So59BtV/biu61PVGcAMgOTkZE1JSfH6e10uFzUpfyHeWb2f5775lt7tm/PqDwbTtnljn9bnDX+0OxA51e6We1sCOPae2+cdWnzZbp8lDlWt7b2HmUD56zjxnn31iqryzMLtPLdkJ5d3j+HF2wYSZTPZhqT/N+b/OR2CMXUikH+DrQUSRaQL7oRxM3CrsyHVjKryp0+38cqKPdyYHM+fJvUlItyRG9lMALg04VKnQzCmTjh1O+4kEckAhgGfisgCz/72IjIPQFVLgGnAAmAbMFtVtzgR74VQVf7wyVZeWbGHKZd25i/X97OkEeJWpq9kZfpKp8MwptacuqvqI+CjSvZnAePLbc8D5vkxtDqhqjw2ZwtvrNrHXcO78OiEXvZ8huHXi38N2Hocpv4L5EtV9dYzC7fzxqp93DOiC78eb0nDGBNc7NpJHXt3zX6eW7KTm5ITLGkYY4KSJY46tCT1EI/8ZzMpPWJ4fFIfSxrGmKBkiaOObMk6xv1vf0PvuOa8cOtAGwg3xgQtG+OoA8cLT3P/21/TokkEr04ZbM9pmEo9O+5Zp0Mwpk7Yb7haUlUe/mAT6UdPMWvqUGKaNXI6JBOgbDp1EyzsekotvbFyL/O+Pcgvx/UguXNrp8MxAWzR7kUs2r3I6TCMqTXrcdTCpow8/jRvG1f0ass9I7o6HY4JcI8vfxywlQBN/Wc9jgt0urSMX7y/idZRDXn6hiS7g8oYEzKsx3GBXvtyD6kHT/DP2wfRIjLC6XCMMcZvrMdxAdKPnOSZhTu4olcs37s41ulwjDHGryxx1JCq8tuPNyMCv7/2YrtEZYwJOXapqobmbz7I0rQcHrm6Fx1aNnE6HFOPvDzhZadDMKZOWOKogdIy5a+fp9E9tilTLu3sdDimnukR3cPpEIypE3apqgY+23yAXTkF/GRMIg1sShFTQ3PT5jI3ba7TYRhTa9bj8FJZmfKPxTu5KCaKq/rEOR2OqYeeXvU0ANf0uMbhSIypHadWALxBRLaISJmIJJ+n3F4R+VZENojIOn/GWNHnWw+RdugE/zs6kfAwGxA3xoQup3ocm4HrAG9GC0epaq6P4zkvVeUfS3bQuU0kE/pZb8MYE9oc6XGo6jZVTXOi7guxJDWbLVnHuX9UNxvbMMaEPFFV5yoXcQEPqWqll6FEZA9wFFDgZVWdcZ5zTQWmAsTGxg6aOXOm13Hk5+fTtGnTKo//Zc0psk8qfxnZhAZBdJmqunYHK6faPX3DdACeTXJmenX7vENLTds9atSo9apa5dDBd6iqT76ARbgvSVX8urZcGReQfJ5zdPD82xbYCIz0pu5BgwZpTSxdurTKYwfyTmnnhz/RZxam1eic9cH52h3MnGr3/rz9uj9vvyN1q9rnHWpq2m5gnXr5+91nYxyqWuspQFU10/Nvtoh8BAwBltf2vDXxyaYsVGFi//b+rNYEoYQWCU6HYEydCNgL9iISJSLNzrwGrsTdY/Grjzdk0S++BV1jQq+ra+rWrM2zmLV5ltNhGFNrTt2OO0lEMoBhwKcissCzv72IzPMUiwVWiMhGYA3wqarO92ecu3Ly+TbzmPU2TJ14ad1LvLTuJafDMKbWHLkdV1U/Aj6qZH8WMN7zejfQ38+hfcecDVmIwDWWOIwx5qyAvVTlNFVlzsYshnVtQ2zzxk6HY4wxAcMSRxW+zTzGntwCrk2y3oYxxpRniaMKH2/IomF4GOMutifFjTGmPJvksAqLth1iRGK0LQtr6sz7N77vdAjG1AnrcVTiaEEx+w6fJLlza6dDMUEkOjKa6Mhop8MwptYscVRiU+YxAPrHt3A4EhNMXt/wOq9veN3pMIypNUscldiYngdAH0scpg5Z4jDBwhJHJTZl5NE1JormjW18wxhjKrLEUYGqsjHjGEnxLZ0OxRhjApIljgoOHi8k50QR/ewylTHGVMoSRwUb090D4/0SrMdhjDGVsec4KtiYkUeDMKF3XHOnQzFBZt5t86ovZEw9YImjgk0ZefRo14zGEeFOh2KCTGREpNMhGFMn7FJVOWVlyqaMY/S3y1TGB15c+yIvrn3R6TCMqTVLHOXsPVzAicISe/DP+MTsLbOZvWW202EYU2uWOMrZlOEZGLdbcY0xpkpOrQD4lIikisgmEflIRCr9TS0i40QkTUR2isjDvo5rY0YejSPCSGxry8QaY0xVnOpxLAT6qGo/YDvwq4oFRCQceAG4CugN3CIivX0Z1Mb0PPp2aEGDcOuIGWNMVRz5Damqn6tqiWfzKyC+kmJDgJ2qultVi4GZwLW+iul0aRlbso7bZSpjjKmGqKqzAYjMBWap6r8r7J8MjFPVH3m27wAuUdVpVZxnKjAVIDY2dtDMmTO9jiE/P5/IqCjST5TROFyIjQqNHkd+fj5Nm4beZTlrd2ixdntn1KhR61U12ZuyPnuOQ0QWAe0qOfQbVf3YU+Y3QAnwdm3rU9UZwAyA5ORkTUlJ8fp7XS4XNSkfLKzdocXaHVp82W6fJQ5VveJ8x0VkCjABGKOVd3sygYRy2/GefcYYYxzk1F1V44BfABNV9WQVxdYCiSLSRUQaAjcDc/wVozHGmMo5dTH/eaAZsFBENojIPwFEpL2IzAPwDJ5PAxYA24DZqrrFoXiNMcZ4ODJXlap2q2J/FjC+3PY8wGaGM8aYABIatw8ZY4ypM5Y4jDHG1IglDmOMMTViicMYY0yNOP7kuC+ISA6wrwbfEg3k+iicQGbtDi3W7tBS03Z3UtUYbwoGZeKoKRFZ5+2j9sHE2h1arN2hxZfttktVxhhjasQShzHGmBqxxOE2w+kAHGLtDi3W7tDis3bbGIcxxpgasR6HMcaYGrHEYYwxpkZCJnGIyDgRSRORnSLycCXHG4nILM/x1SLS2f9R1j0v2v2AiGwVkU0islhEOjkRpy9U1/Zy5a4XERWRoLhl05t2i8iNns99i4i84+8YfcGLn/WOIrJURL7x/LyPr+w89YmIvCoi2SKyuYrjIiLPed6TTSIysE4qVtWg/wLCgV1AV6AhsBHoXaHMfcA/Pa9vxr2creOx+6Hdo4BIz+t7g6Hd3rbdU64ZsBz4Ckh2Om4/feaJwDdAK892W6fj9lO7ZwD3el73BvY6HXcdtHskMBDYXMXx8cBngABDgdV1UW+o9DiGADtVdbeqFgMzgWsrlLkWeMPz+n1gjIiIH2P0hWrbrapL9b+LaX2Fe6XFYODNZw7wR+AvQKE/g/Mhb9p9D/CCqh4FUNVsP8foC960W4HmntctgCw/xucTqrocOHKeItcCb6rbV0BLEYmrbb2hkjg6AOnltjM8+yoto+5FpI4BbfwSne940+7y7sb910kwqLbtnm57gqp+6s/AfMybz7w70F1EvhSRrzwrctZ33rT7MeB2EcnAvc7P//onNEfV9HeAVxxZyMkEHhG5HUgGLnc6Fn8QkTDgb8AUh0NxQgPcl6tScPcwl4tIX1XNczQq37sFeF1VnxaRYcBbItJHVcucDqy+CZUeRyaQUG473rOv0jIi0gB3V/awX6LzHW/ajYhcAfwG9xrwRX6Kzdeqa3szoA/gEpG9uK//zgmCAXJvPvMMYI6qnlbVPcB23ImkPvOm3XcDswFUdRXQGPdEgMHMq98BNRUqiWMtkCgiXUSkIe7B7zkVyswBfuB5PRlYop7RpXqs2naLyADgZdxJIxiudZ9x3rar6jFVjVbVzqraGff4zkRVXedMuHXGm5/1/+DubSAi0bgvXe32Z5A+4E279wNjAESkF+7EkePXKP1vDnCn5+6qocAxVT1Q25OGxKUqVS0RkWnAAtx3X7yqqltE5A/AOlWdA/wf7q7rTtyDTTc7F3Hd8LLdTwFNgfc89wLsV9WJjgVdR7xse9Dxst0LgCtFZCtQCvxcVet179rLdj8I/EtEfoZ7oHxKff/jUETexf1HQLRn7OZ3QASAqv4T91jOeGAncBL4YZ3UW8/fN2OMMX4WKpeqjDHG1BFLHMYYY2rEEocxxpgascRhjDGmRixxGGOMqRFLHMYYY2rEEocxxpgascRhjB+ISGcRSRWRt0Vkm4i8LyKRTsdlzIWwxGGM//QAXlTVXsBx3GvAGFPvWOIwxn/SVfVLz+t/A5c5GYwxF8oShzH+U3F+H5vvx9RLljiM8Z+OnnUgAG4FVjgZjDEXyhKHMf6TBtwvItuAVsBLDsdjzAUJiWnVjQkQJap6u9NBGFNb1uMwxhhTI7YehzHGmBqxHocxxpgascRhjDGmRixxGGOMqRFLHMYYY2rEEocxxpga+f9r8QMZ5Ofj7gAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(p, np.log10(odds),label=\"log10odds\")\n",
    "plt.vlines(0.5, -2, 2, color=\"g\", linestyle=\"dashed\")\n",
    "plt.hlines(0, 0.4, 0.6, color=\"g\", linestyle=\"dashed\")\n",
    "plt.text(0.5, 0, \"(0.5, 0)\")\n",
    "plt.xlabel(\"p\")\n",
    "plt.ylabel(\"log10odds\")\n",
    "plt.legend()\n",
    "plt.title(\"log10odds\")\n",
    "plt.grid()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "对数几率的符号也可以划分类别，分类器可以为$C=sign(logit)$"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python [conda env:anaconda3]",
   "language": "python",
   "name": "conda-env-anaconda3-py"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.8"
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": true,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": false,
   "toc_position": {},
   "toc_section_display": true,
   "toc_window_display": true
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
